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is  sought  through  fragment  analyses,  decompositions  of  the  density  of  states  and  o\erl; 
population  studies.  All  fragments  tend  to  restore  their  missing  C-H  bonds  when 
bound  to  these  surfaces  -  CHI’  prefers  the  on-top,  CH  '  the  bridging  and  CH  the 
capping  geometry.  CH^  anchors  more  strongly  to  the  on-top  site  of  a  metal  surface 
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BONDING  AND  COUPLING  OF  C  FRAGMENTS  ON  METAL  SURFACES. 

Chong  Zheng  ' ,  Yitzhak  Apeloig  ,  and  Roald  Hoffmann*4 

Contribution  for  athe  Department  of  Chemistry  and  Materials  Science 
Center,  Cornell  University,  Ithaca,  NY  14853;  ^the  Department  of 
Chemistry,  Technion-Israel  Institute  of  Technology,  Haifa,  30012  Israel. 

Abstract;  The  bonding  and  reactivity  of  CH^,  CH2  and  CH  fragments  to  Ti 
(0001),  Cr(110)  and  Co(0001)  metal  surfaces  is  examined  with  extended  Kuckei 
band  calculations  on  two  dimensional  slabs  of  metal  and  adsorbate.  A  local 
chemical  viewpoint  is  sought  through  fragment  analyses,  decompositions  rf 
the  density  of  states  and  overlap  population  studies.  All  fragments  tend 
restore  their  missing  C-H  bonds  when  bound  to  these  surfaces  -  CH  prefers 
the  on-top,  CH^  the  bridging  and  CH  the  capping  geometry.  CH^  anchors  more 
strongly  to  the  on-top  site  of  a  metal  surface  of  high  d  band  filling  since 
the  antibonding  feature  at  the  top  of  the  d  band  that  results  from  j-bonding 
destabilizes  sites  of  higher  coordination.  Similar  conclusions  holu*  for 
other  fragments.  Thus,  the  mobility  of  these  fragments  is 
reduced  on  metal  surfaces  of  higher  d  band  filling.  The  mobility  patterns 
of  CH^,  CH^ ,  CH  are  examined.  In  general,  on  the  way  to  products  there  are 
barriers  to  migration  on  the  surface,  a  proximity  or  crowding  effect  which 
makes  it  costly  for  two  fragments  to  approach  on  the  surface,  and  a  barrier, 
small  or  large,  to  their  reaction  with  each  other.  When  two  Cx  fragments 
couple, the  C-C  a*  orbital  rises  from  below  the  Fermi  level.  It  is  initially 
filled,  then  empties  as  the  reaction  proceeds.  Hence  the  lower  the  Fermi 
level  (for  metals  at  the  right  side  of  the  transition  series),  the  smaller 
the  reaction  barrier.  The  decrease  of  the  mobility  and  the  lower  coupling 
barrier  as  the  metal  is  changed  from  the  left  to  the  right  side  in  the  Period 
Table  may  be  two  of  many  reasons  why  metals  in  the  middle  of  the  transition 
series  have  higher  reactivity  in  Fischer-Tropsch  catalysis. 

Present  address;  Department  of  Chemistry,  Stanford  University, 

Stanford,  CA  94305 
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The  Fischer-Tropsch  (FT)  synthesis,  which  cen  be  defined  ee  the  reductive 

oligomerization  of  carbon  monoxide  over  a  heterogeneous  catalyst  (eq.  1),  was 

described  nearly  60  years  ago.  '  Because  of  the  great  technological  importance  o: 

this  reaction  much  effort,  especially  during  the  last  two  decades,  has  been  devoted 

to  the  elucidation  of  its  mechanism. Although  the  subject  still  continues  to  :e 

Cotolyst 

n  CO  ♦  mH,  - —  C,  H,  Of  (I) 

strongly  debated4*  the  accumulated  evidence  suggests  that  under  conditions  which  lest 
to  oxygen- free  products  (l.e.,  z  -  0  in  eq.  1  the  FT  reaction  proceeds  via  the 
"carbide/methylene"  mechanism^* ^  which  is  drawn  schematically  in  Scheme  1.  Under 
these  conditions,  the  major  products  era  a-olefins  and  hydrocarbons,  and  the  oxyge- 
ends  up  primarily  as  water,  along  with  soma  oxygenated  products. 
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Scheme  1 

The  carblda/carbana  mechanism  was  first  suggested  by  Fisher  and  Tropsch 
themselves  as  early  as  1926,*  and  it  was  re- Introduced  with  additional  details  by 
Craxford  and  Rldeal.3  According  to  this  mechanism  carbon  monoxide  is  first  adsorbed 
and  then  dissociates  on  the  metal  surface  to  give  "surface  carbides”  (step  1,  Scheme 
1),  which  are  then  hydrogenated  to  give  surface -bonded  methylene  and  methyl 
fragments  (step  2).  The  oligomerization  of  these  metal-bonded  fragments  (step  3)  is 
followed  by  a  termination  step  such  as  a  ^-elimination  (step  A),  which  followed  bv 
desorption  yields  the  final  products.^ 

The  most  extensive  and  convincing  support  for  the  carbide/methylene  mechanism 
comes  from  the  elegant  studies  by  the  research  groups  of  Biloen5  and  Petit6. 
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The  Fischer-Tropsch  (FT)  synthesis,  which  can  be  defined  as  the  reductive 
oligomerization  of  carbon  monoxide  over  a  heterogeneous  catalyst  (eq.  1),  was 
described  nearly  60  years  ago.  *  Because  of  the  great  technological  importance 
this  reaction  much  effort,  especially  during  the  last  two  decades,  has  been  devoted 
to  the  elucidation  of  its  mechanism.1"3.  Although  the  subject  still  continues  to  be 


:e  or 


„  Catalyst 
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(I) 


strongly  debated  the  accumulated  evidence  suggests  that  under  conditions  which  lead 
to  oxygen- free  products  (i.e.,  z  ■  0  in  eq.  1  the  FT  reaction  proceeds  via  the 
"carbide/methylene"  mechanism1"3  which  is  drawn  schematically  in  Scheme  1.  Under 
these  conditions,  the  major  products  are  a-olefins  and  hydrocarbons,  and  the  oxvger, 
ends  up  primarily  as  water,  along  with  some  oxygenated  products. 
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The  carbide/carbene  mechanism  was  first  suggested  by  Fisher  and  Tropsch 
themselves  as  early  as  1926,  and  it  was  re-introduced  with  additional  details  bv 
Craxford  and  Rideal.3  According  to  this  mechanism  carbon  monoxide  is  first  adsorbed 
and  then  dissociates  on  the  metal  surface  to  give  "surface  carbides'*  (step  1,  Scheme 
1),  which  are  then  hydrogenated  to  give  surface -bonded  methylene  and  methyl 
fragments  (step  2).  The  oligomerization  of  these  metal-bonded  fragments  (step  3)  is 
followed  by  a  termination  step  such  as  a  ^-elimination  (step  4),  which  followed  bv 
desorption  yields  the  final  products.1 

The  most  extensive  and  convincing  support  for  the  carbide/methylene  mechanism 
comes  from  the  elegant  studies  by  the  research  groups  of  Biloen3  and  Petit^  Biloer. 
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and  Sachtler  found  that  Ni,  Co  and  Ru-based  catalysts  which  are  pre -doped  with  'L'Z 


labeled  carbon  yield  upon  treatment  with  JC0/H2  mixtures  under  FT  conditions  a 


product  mixture  consisting  mainly  of  JCH^  and  of  hydrocarbons  containing  several 
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C-atoms  in  the  same  molecule  .  These  results  indicate  that  the  carbidic  species 


once  formed,  can  react  with  hydrogen  to  give  CHX  intermediates  (x  -  1-3)  which 


polymerize  to  produce  hydrocarbons,  in  agreement  with  steps  2  and  3  in  Scheme  1. 


The  feasibility  of  the  first  step  in  the  mechanism,  in  which  the  carbidic  surface 


is  formed,  has  been  demonstrated  convincingly  by  other  investigators.7  Thus,  a 


rapid  dissociative  chemisorption  of  carbon  monoxide  has  been  shown  to  occur 


various  clean  metal  surfaces,  including  typical  FT  catalysts  such  as  iron.7  More 


importantly,  the  formation  of  a  "carbidic  layer"  under  real  FT  conditions  is  well 


known. ' 


The  beautiful  studies  of  Petit  and  Brady  on  the  induced  decomposition  of 


gaseous  diazomethane  (CHjNj)  on  typical  FT  catalysts  provides  additional  strong  ar.d 


independent  evidence  for  the  operation  of  the  carbide/carbene  mechanism.1 


Decomposition  of  on  Ni- ,  Pd-,  Fe-,  Co-,  Ru- ,  and  Cu-  surfaces,  at  atmospheric 


pressure  and  in  the  temperature  range  of  25-250  °C,  produces  only  ethylene  and 


dinitrogen.  This  indicates  that  in  the  absence  of  hydrogen  the  absorbed  CH- 


fragments  dimerize  to  ethylene  (Scheme  2) ,  but  polymerization  to  higher  hydrocarbons 


CH,  CH. 

/\_  n* 


CHj-CH* 

T77777T77T 


CHj-CH, 
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Scheme  2 


ses  not  occur.  However,  reaction  of  a  mixture  of  Hj  and  over  Co-,  Fe-,  and 


Ru- surfaces,  all  typical  FT  catalysts,  produces  a  variety  of  hydrocarbons  with 


isomer  and  molecular  weight  distribution  typical  of  a  "real"  FT  reaction. 


Furthermore,  decomposition  of  CH2N2  on  surfaces  which  are  not  capable  of 
dissociative  chemisorption  of  H2,  such  as  Cu,  yields  only  ethylene,  even  in  the 


presence  of  H?.  Wang  and  Ekerdt  in  a  more  recent  study  showed  that  pyridine  can 


used  to  scavenge  C^-Cj  alkyl  species  from  the  surface  of  an  iron  catalyst  during  : 

Q 

synthesis.  These  findings  are  also  consistent  with  the  carbide/carbene  mechanise 

in  which  alkyl  fragments  are  the  immediate  precursors  to  the  FT  products.  In 

another  important  recent  paper,  Ekstrom  and  Lapszewicz  showed  that  high  molecular 

weight  hydrocarbons  can  be  formed  by  the  reaction  of  carbides  with  hydrogen  in  the 

presence  of  water^.  Furthermore,  Winograd  et  al  have  recently  reported  direct 

observation  of  CH,  CH2  and  CH^  intermediates  on  a  Ni(lll)  methanation  catalyst.  “ 

Additional  support,  although  indirect,  for  the  chemical  feasibility  of  the 

various  steps  of  the  carbide/methylene  mechanism  stems  from  the  study  of  appropria 

model  organometallic  complexes.  Substantial  experimental  effort  has  been  devoted 

la 

the  last  decade  to  the  study  of  such  model  complexes  Relevant  systems  and 
molecular  reactions  will  be  mentioned  briefly  along  the  paper,  and  reviews  (in 
addition  to  ref.  la)  are  gathered  in  reference  11. 

Despite  the  extensive  study  of  the  FT  synthesis  many  of  the  mechanistic  detai 
remain  poorly  understood.  In  particular,  there  is  virtually  no  knowledge  of  the 
electronic  and  geometrical  factors  that  control  the  polymerization  process  (step  2 
Scheme  1)  and  very  little  is  known  of  the  nature  of  the  bonding  of  the  organic 
fragments  to  the  surface.  Representative  questions  that  remain  unanswered  are: 
What  are  the  most  effective  binding  sites  on  the  metal  surface  for  the  adsorbed 
fragments  and  are  these  sites  the  same  for  all  radicals?  What  is  the  favored 
orientation  for  the  coupling  of  two  methylenes  to  produce  ethylene  (Scheme  2)?  Ar 
the  binding  sites  and  the  dimerization  mechanism  the  same  for  different  metals? 
Vhat  is  the  relative  mobility  of  CH2  and  CH3  fragments  on  the  metal  surface,  and  h 
does  this  mobility  change  as  the  metal  is  varied?  At  present  these  and  similar 
questions  are  difficult  to  answer  experimentally,  although  encouraging  advances  ha 
been  made  recently. ^  We  will  address  these  and  other  interesting  mechanistic 
questions  regarding  the  FT  synthesis  in  this  paper. 

A  full  theoretical  treatment  of  the  FT  reaction  is  a  vast  project  and  we  have 
to  impose  some  limitations  on  the  scope  of  problems  that  will  be  tackled.  In  this 
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study  we  concentrate  on  three  major  issues:  (1)  The  mode  of  binding  of  the 


postulated  reaction  intermediates  --  methyl,  methylene  and  methyne  --  to  the  metal 


surface.  (2)  The  mobility  of  these  radicals  on  the  metal  surface.  (3)  The 


electronic  and  geometrical  requirements  for  the  surface  induced  coupling  of  the 


organic  radicals. 


We  use  in  our  work  tight  binding  extended  Huckel  calculations,^  with  details 


given  in  Appendix  I.  The  analysis  of  the  surface  calculations  is  based  on  the 


methodology  and  the  tools  that  we  have  described  in  detail  in  our  recent  study  of 


H  and  H-H  activation  on  surfaces.  3  The  extended  Hiickel  method,  whether  applied  to 


discrete  molecules  or  extended  systems,  has  well-known  limitations.  It  does  not  do 


well  at  predicting  distances,  and  this  will  impose  severe  limitations  in  the  analses 


of  mobility  and  coupling.  But  the  method  does  seem  to  capture  the  essence  of 


bonding.  It  is  also  transparent  and  useful  in  constructing  explanations,  and  it  is 


for  this  reason  that  we  use  it. 


While  many  of  the  problems  we  will  attack  have  not  been  investigated 


:heoretically  before,  some  have.  The  relevant  studies  of  Baetzold,  Muetterties 


Shustorovich , ^ ^  of  Minot,  Van  Hove  and  Somorjai,^®  of  Anderson, as  well  as 


other  related  theoretical  work  will  be  discussed  in  the  paper.  We  also  note  that 


the  first  step  of  the  FT  reaction,  the  chemisorption  of  CO  on  various  metal 


surfaces,  has  been  studied  extensively  by  many  theoreticians . ^ • 18b-d, 19 


This  will  be  a  long  paper.  It  could  easily  have  been  chopped  up  into  three 


papers,  one  on  fragment  bonding,  a  second  on  migration,  a  third  on  fragment 


coupling.  We  think  nothing  would  be  gained  by  this,  and  continuity  lost.  We  ask 


the  reader  to  bear  with  us  as  we  analyze  In  detail,  an  important  reaction  at  the 


border  of  chemistry,  physics  and  catalysis. 


S  V 


The  FT  synthesis  occurs  on  Fe  and  Co  catalysts  and  to  some  extent  also  on  Ni 

surfaces.^-  A  major  goal  of  any  mechanistic  study  of  the  FT  synthesis  is  to  cry  to 

define  and  understand  the  factors  that  determine  the  catalytic  reactivity  of  the 

metal  surface.  In  this  respect  two  key  parameters  are  the  identity  of  the  metal  and 

the  lattice  form  in  which  it  crystallizes . ^ ^  For  each  metal,  one  can  choose 

surfaces  with  different  indices,  and  generally  each  surface  is  expected  to  exhibit  a 

different  reactivity.  We  are  going  to  use  a  fixed  surface  throughout  the  paper,  to 

model  what  must  be  only  one  aspect  of  reality.  The  surface  we  chose  is  the  (0001) 

surface  of  a  hexagonal  metal.  With  the  high  symmetry  of  this  surface  the 

computational  times  can  be  greatly  reduced.  For  example,  when  a  methyl  group  is 

adsorbed  on  this  surface  in  an  "on- top"  or  "three -fold"  geometry,  the  hexagonal 

symmetry  is  still  retained  and  the  computational  advantage  of  this  over  other 

surfaces  of  FT  catalysts  (e.g.  Fe  (110))  is  obvious.  The  results  can  also  be 

15  71 

compared  with  our  previous  studies  on  other  hexagonal  surfaces.  ’ 

To  extract  the  basic  electronic  effects  that  determine  the  reactivity  of  the 

metal,  we  are  going  to  compare  three  metal  surfaces:  Ti  (0001),  Cr  (110)  and  Co 

(0001).  They  are  shown  in  1,  2  and  3,  respectively.  Although  the  Cr  (110) 

surface  is  not  hexagonal,  it  can  be  viewed  as  a  distortion  therefrom. 

Ti  (0000-  Cr  (1 10)  Co  (OOOI) 


H  1  h  H*.si 

I  2  3 


Ley  us  begin  with  the  hep  lattice  of  cobalt  and  employ  the  experimental  Co-Co 
distance  of  2.51  A.  The  next  choice  to  be  made  is  that  of  the  thickness  of  the  slab 
of  the  metal  to  be  used  in  the  calculations,  of  necessity  a  compromise  between 
computational  economy  and  reasonable  accuracy.  Our  previous  studies^0  ’  ^a’  ^ 

as  well  as  exploratory  calculations  on  the  the  cobalt  slab  lead  to  the  conclusion 


that  a  reasonable  choice  is  a  3 -layer  slab.  We  use  here  a  slab  of  three  layers 
because  the  changes  in  the  important  surface  properties  (e.g.,  the  Fermi  level, 
charge  distribution,  overlap  population  etc.)  are  small  on  going  to  a  4-layer  slab. 
The  3 -layer  slab  model  for  the  hep  Co  surface  is  shown  in  a  side  view  in  4.  We 
choose  the  z  axis  to  be  perpendicular  to  the  surface.  A  top  view  is  in  5 ,  where 
the  dots  in  the  triangular  hollows  represent  the  Co  atoms  in  the  next  layer  below. 


4,  tide  S,  top 

Figure  1  shows  the  density  of  states  (DOS)  curve  of  the  slab.  A  DOS  curve  is 
the  solid  state  analogue  of  an  energy  level  diagram  and  it  gives  the  number  of 
levels  in  a  particular  energy  interval.  The  metal  bands  (or  orbitals)  are  filled  u 
to  the  Fermi  level  --  the  solid  state  equivalent  of  the  molecular  HOMO.  The  dashed 
lines  in  Figure  1  refer  to  the  total  DOS  curve.  The  darkened  areas  indicate  the 
contribution  of  the  s  (Figure  la)  and  the  p  states  (Figure  lb)  to  the  total  DOS. 

The  darkened  areas  are  examples  of  projected  or  local  DOS  curves  which  single 


out  the  contribution  of  a  certain  atom  or  a  group  of  atomic  or  fragment  orbitals  to 
the  overall  DOS  plot.  The  states  that  are  not  s  or  p  are  d  states.  The  dotted 
lines  are  integration  curves,  from  0  to  100%,  which  additively  count  the  relative 
number  of  states  occupied  as  one  sweeps  up  the  energy  scale.  It  is  clear  from 
Figure  1  that  a  substantial  number  of  s  and  some  p  states  penetrate  into  the  d  band 
On  the  average  any  Co  atom  has  its  s  band  approximately  one  third  filled. 

In  any  reactivity  problem,  molecular  or  solid  state,  the  energy  and  bonding 
capabilities  of  the  frontier  orbitals,  the  lower  unoccupied  and  higher  occupied 
orbitals  of  the  system  will  play  a  crucial  role.  In  the  case  of  a  metal  surface 
these  are  the  orbitals  near  the  Fermi  level.  It  is  important  to  know  how  the 
essential  features  of  these  frontier  orbitals  of  the  surface  change  as  the  metal  is 


1  wi 


varied.  As  the  number  of  the  d  electrons  increases  from  Ti  to  Co,  the  nuclear 
charge  is  less  effectively  "screened"  and  d  electrons  are  more  strongly  bound  by  the 
bigger  Coulomb  interaction.  There  are  therefore  two  factors  competing  in  the 
determination  of  the  Fermi  level:  the  filling  of  the  d  band,  which  tends  to  raise 
the  Fermi  level,  and  the  Coulomb  interaction  that  pulls  down  the  Fermi  level  with 
increasing  d  electron  count.  It  turns  out  that  the  Coulomb  interaction  wins  out, 
and  the  Fermi  level  descends  slowly  along  the  right  side  of  the  transition  series 
(the  relevant  work  functions  are  Ti:  4.33,  Cr:  4.5,  Co:  5.0  eV^®d) .  Another 
consequence  of  the  increasing  Coulomb  attraction  is  the  decreasing  d  band  width 
along  the  transition  series,  due  to  the  "tighter"  wave  functions  of  the  d  electrons. 
The  calculations  confirm  these  general  considerations.  The  extended  Huckel  method, 
as  usaual ,  gives  a  much  too  high  magnitude  of  the  ionization  potential  or  Fermi 

level.  It  does  so  for  molecule  as  well  as  for  extended  structures. 

IS  21 

In  analogy  to  our  previous  studiesiJ’  .  we  find  that  the  charge  distribution 
among  the  slab  layers  changes  as  a  function  of  the  number  of  electrons  per  metal 
atom.  Basically,  this  results  from  the  fact  that  an  inner  atom  has  more  neighbors 
than  a  surface  one.  The  band  dispersion  is  a  function  of  inter-unit-cell 
interaction,  so  that  the  more  interactions  (neighbors)  one  has,  the  wider  the 
resulting  band.  Thus  the  states  of  the  "surface  atoms"  (layer  A  in  4)  form 
narrower  bands  than  the  bulk- like  atoms  (layer  B  in  4).  The  wider  bulk  bands  are 
filled  first  and  thus  the  bulk  atoms  become  negative  relative  to  the  surface.  At 
some  point  along  the  transition  series  the  two  layers  will  have  equal  charges  and 
past  this  point  the  surface  layer  will  become  negative.  This  argument  is  presented 
schematically  in  Figure  2.  The  conclusions  from  these  qualitative  arguments  are 
that  at  the  right  side  of  the  transition  series  surfaces  are  expected  to  be  negative 
relative  to  the  bulk  (or  to  the  isolated  atom),  while  at  the  left  side  of  the 
transition  series  surface  are  expected  to  be  positive. 
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Figure  2:  A  schematic  picture  shoving  the  relative  charge  on  surface  and  bulk  atoms 
(b)  caused  by  different  effective  band  width  for  each  kind  of  atoms  (a) . 
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The  calculations  fully  support  these  conclusions.  For  Co  we  find  that  each  of 

the  surface  atoms  carries  a  negative  charge  of  0.305  electrons  (note  that  because  : : 

the  ABA  arrangement  the  charge  on  each  bulk  metal  is  twice  as  large  and  of  opposite 

sign  to  charge  on  a  surface  atom).  For  Cr  we  get  a  negative  charge  of  -0.288  at  a 

surface  atom.  For  Ti  the  polarization  is  small  and  both  surface  and  bulk  layers  are 

nearly  neutral  (each  surface  Ti  atom  carries  a  charge  of  -0.025  electrons).  To  put 

these  charge  distributions  in  a  chemical  context  we  might  say  that  the  surface 

layers  of  the  Co  and  Cr  slabs  can  be  described  as  being  nucleophilic  and  Ti  neutral. 

Better  calculations  will  temper  the  indubitably  exaggerated  density  shifts  between 

surface  and  bulk,  but  the  trend  should  remain. 

Another  important  difference  between  the  surface  and  the  bulk  atoms,  shown  in 

Table  1,  is  in  the  electron  distribution  between  the  atomic  metal  orbitals.  Let  us 

examine  first  the  surface.  The  data  in  Table  1  reveal  in  the  Co  case,  a  significant 

electron  flow  from  the  "in-plane"  orbitals  (dv2  v2,  d„„)  into  orbitals  which  are 

x.  •y  xy 

perpendicular  to  the  surface  (dz2,  dxz ,  dyz)  in  the  Co  case.  But  as  the  d  electron 
filling  decreases  the  situation  is  reversed:  for  Ti  there  is  an  electron  flow  from 
dz2,  dxz  and  dyz  to  dx2-y2  and  dxy  T*11-3  again  due  to  the  fact  the  "in-plane" 
orbitals  on  the  surface  overlap  better  with  their  neighbors,  resulting  a  wider  band. 
At  low  electron  counts  these  orbitals  are  filled  first,  but  as  the  filling  increases 
the  narrower  bands,  which  are  at  intermediate  energy,  are  filled  more. 

Table  _i 

For  the  bulk  orbitals,  the  perpendicular  orbitals  (dz2,  dxz.dyZ)  overlap  not 
only  with  orbitals  in  the  same  plane  (B  in  4)  but  also  the  layers  sandwiching  them 
(A  layers  in  4),  thus  the  band  widths  are  bigger  than  those  of  d  2  2  and  dv„. 

~  x  -y  xy 

For  the  same  reason  that  we  have  adduced  above  the  charge  flows  from  dz2,  dxz  and 
d  to  d  2  2  and  d  at  high  electron  counts,  but  is  reversed  for  low  electron 

yz  x  *y  xy 

fillings . 

We  are  now  in  a  position  to  bring  to  the  surface  a  layer  of  organic  molecules 


TABLE  1.  Electron  Distribution  Among  the  Metal  Atomic  Orbitals  of  Co (0001) , 
Cr(110),  and  Ti(0001)  Slabs.  Coordinate  Axes  are  given  in  1. 


Co 

d9 

Cr 

d6 

Ti 

d4 

Orbital 

"bulk" 

"surface" 

"bulk" 

"surface" 

"bulk" 

d  2  2  or  d 

x  -y  xy 

1.569 

1.541 

1.034 

0.896 

0.632 

0.566 

V 

1.731 

1.510 

1.041 

0.947 

0.566 

0.570 

d  or  d 

xz  yz 

1.764 

1.445 

1.046 

0.856 

0.588 

0.612 

s 

0.651 

0.616 

0.752 

0.666 

0.798 

0.730 

P 

0.258 

0.293 

0.334 

0.305 

0.221 

0.295 

Total 

9.305 

8.390 

6.288 

5.424 

4.025 

3.951 
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METAL  SURFACES 


In  order  to  study  the  adsorption  of  CHj  (x  -  1-3)  radicals  on  the  metal  surras 
we  have  to  simplify  the  calculations  further  by  covering  only  one  side  of  the  2- 
iayer  slab.  For  justification  we  rely  on  our  previous  study  which  showed  thac  verv 
similar  results  are  obtained  for  a  coverage  of  a  4- layer  Ni  slab  with  H2  on  one  si cL 
or  on  both  sides. ^ 

The  one-side  coverage  makes  the  two  identical  A  layers  in  4  different  in  the 
covered  metal  slab,  as  shown  schematically  in  6.  The  top  layer  A'  can  be 
described  as  the  "adsorbing"  layer,  the  inner  layer  B'  is  "bulk- like"  and  the  bottor 
layer  C'  is  similar  in  character  to  the  surface  layer  A  in  the  bare  metal  (except 
for  a  small  perturbation  by  the  remote  adsorbent). 


We  have  studied  the  adsorption  of  three  organic  radicals,  CHj ,  CH2  and  CH,  all 
believed  to  be  intermediates  in  the  FT  synthesis.^  Regardless  of  the  geometry  that 
is  chosen,  a  1:1  coverage  of  the  metal  surface  by  CH^  or  CH2  is  chemically 
unrealistic  due  to  the  very  short  distances  and  the  resulting  excessive  steric 
repulsions  between  the  hydrocarbons.  We  have  chosen  a  one-third  coverage,  which 
ensures  minimal  interactions  between  neighboring  fragments  but  which  still  allows 
the  use  in  the  calculation  of  a  convenient  unit  cell  which  is  only  three  times 
larger  than  in  the  bare  metal.  The  details  concerning  the  unit  cell  which  was  used 
the  Brill ouin  Zone,  the  special  k  points,  etc.  are  given  in  Appendix  I.  For  each  0 
the  organic  radicals  we  have  considered  three  possible  adsorption  sites:  a  mono  - 
coordinated  "on- top"  site  7,  a  site  "bridging"  two  metal  atoms  8,  and  "triply - 
bridging",  "capping"  or  "hollow"  site  9. 
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As  in  our  previous  studies  *  ’  we  use  the  language  and  formalism  of  simple 

perturbation  theory.  Within  this  framework  the  interaction  of  two  levels, AE,  is 
given  by  equation  2.  The  magnitude  of  the  matrix  element  Hjj  in  the  numerator  is 


AE 


Kls 


r  •  _  r  0 

Ll 


(2) 


related  to  the  overlap  of  the  relevant  orbitals,  and  the  denominator  tells  us  chat 
the  interaction  is  greater  the  more  two  orbitals  come  into  resonance. 

The  consequence  of  orbital  interactions  between  the  metal  bands  and  the 
orbitals  of  the  organic  fragment  can  be  traced  down  and  analyzed  by  examining  the 
DOS  curves.  Contributions  or  projections  of  specific  orbitals  are  particularly 
helpful.  In  general,  a  strong  shift  in  the  position  of  a  particular  fragment 
orbital  as  it  approaches  the  metal  surface  indicates  strong  interaction,  either 
bonding  or  antibonding.  A  small  shift,  on  the  other  hand,  indicates  little 
interaction. 

Another  important  tool  for  the  analysis  is  the  COOP  (for  Crystal  Orbital 
Overlap  Population)  curve,  which  gives  the  relative  number  of  levels  in  a  given 
energy  interval  weighted  by  the  contribution  that  these  levels  make  to  the  overlap 
population  of  a  specified  bond.  In  other  words  a  COOP  curve  allows  us  to  determine 
if  a  collection  of  energy  levels  contributes  to  bonding  or  antibonding  between  two 
atoms  or  fragments.  We  find  that  both  the  DOS  and  the  COOP  curves,  but  in 
particular  their  combination,  are  very  effective  in  analyzing  the  bonding  properties 
of  metal  surfaces.  The  use  of  these  tools  will  be  demonstrated  throughout  the 
paper.  The  interested  reader  is  referred  to  our  previous  studies140 ’ 15a * 19 ’ 21  for  a 
more  complete  discussion. 

We  proceed  now  to  analyze  in  some  detail  the  bonding  between  each  of  the 
organic  radicals  and  the  metal  surface. 


The  orbitals  of  the  methyl  fragment  are  well  known  and  they  are  shown  in  Figur 
3.  *’e  arbitrarily  choose  to  define  the  methyl  fragment  as  a  radical  so  that  at 

infinite  separation  both  the  metal  and  the  radical  are  neutral.  The  HOMO  of  the 
methyl  radical  is  n-CHj,  a  nonbonding  type  orbital  at  -11.75  eV.  This  orbital, 
which  we  will  call  n,  is  singly  occupied.  Lower  in  energy,  at  -15.56  eV,  are  the 
two  degenerate  x-CH-j  orbitals  and  still  lower  is  the  0-CH3  orbital.  The  LUMO 

|g 

consists  of  the  two  degenerate  x-CHj  orbitals  at  4.99  eV. 


There  are  many  interaction  modes  between  the  methyl  radical  and  the  metal  slab 
but  the  perturbation  expression  helps  one  sort  these  out.  In  the  schematic  diagram 
10  we  have  assigned  number  and  ^  to  the  interactions  of  n,  and  x-CH^ 

with  the  unfilled  band  states  of  the  surface  slab,  and  (Jj),  ^),  to  the 

interaction  with  filled  metal  states.  Whether  these  are  overall  stabilizing  or 
destabilizing  is  a  function  of  the  orbital  filling.  The  magnitude  of  the 
interactions  obviously  depends  critically  on  both  the  separation  in  energy  between 
the  interacting  orbitals  and  the  effectiveness  of  the  overlap  in  question.  10 
also  contains  a  qualitative  summary  of  our  expectations,  including  the  anticipated 
direction  of  charge  transfer  that  is  a  consequence  of  each  interaction. 
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A  word  of  explanation  is  needed  on  these  expectations.  Normally  the  focus 
would  be  on  two-electron  bonding  interactions,  for  these  are  both  stabilizing  and 
effective  in  charge  transfer.  But  in  the  case  at  hand  none  of  the  interactions  of 
this  type,  0 

or  ©is  very  good,  because  the  orbitals  in  question  are  far  from 
resonance.  Instead  interactions^  and  (©  of  the  methyl  radical  orbital  become  most 
important.  And  interaction  (©,  typical  of  energetically  ambiguous  and  difficult  to 
analyze  3 -electron  interactions,  will  be  crucial. 

In  addition  to  these  interactions,  all  of  which  have  molecular  equivalents,  we 
have  interaction©)  characteristic  of  surfaces,  ^^represents  schem-*-’ tally  the 
metal  slab's  ability  to  shift  electron  density  between  the  b  .Ik  and  surf  .  rers, 
or  on  the  surface  to  shift  density  between  those  metal  atoms  involved  in  l  '-'ng  and 
those  left  alone.  This  occurs  in  response  to  the  electronic  demands  imposed  by 
interactions  with  the  organic  adsorbate. 

We  can  trace  the  validity  of  the  perturbation  theory  based  characterization  of 
the  primary  interactions  by  examining  the  consequences  or  symptoms  of  interaction  in 
the  DOS  curves  of  Figure  4a.  This  is  for  methyl  on  Co(0001)  in  the  on-top  geometry. 
Projections  of  the  methyl  orbitals  are  darkened  and  the  position  of  these  orbitals 
in  the  isolated  organic  fragments  are  indicated  by  the  arrows. 


Figure  4:  Total  DOS  (dashed  line)  and  the  CH^  contribution  (darkened  area)  when 
CH3  group  is  chemisorbed  In  an  on-top  geometry  on  Co(0001)  (a)  and  Ti(0001)  (b) 
surfaces.  The  arrows  Indicate  the  CH3  MO  levels  before  the  adsorption  occurs. 


As  expected,  the  o-CH^  and  the  jt-CH^  bands  of  the  adsorbed  CH^  are  essential  1 
at  the  same  energy  as  in  the  isolated  radical  and  their  bands  are  narrow  --  a  clea 
indication  of  their  weak  interaction  with  the  surface.  In  contrast,  the  energies 
the  n  and  the  n  -CH^  orbitals  change  significantly  upon  interaction  with  the  metal 
The  n  band  is  pulled  down  to  a  lower  energy  (by  -  1  eV) ,  while  the  jr*-CHj  band  is 
pushed  up  in  energy  approximately  to  the  same  degree.  The  large  shifts  of  these 
orbitals  relative  to  the  orbitals  of  the  isolated  radical  are  clear  indications  of 
their  strong  interaction  with  the  surface.  The  interaction,  of  course,  depends  or. 
the  M-C  distance.  The  DOS  in  Figure  4  are  calculated  for  a  typical  M-C  distance  o 
2.1  A  (for  the  choice  of  the  M-C  distance  see  Appendix  I). 

Figure  4b  is  the  DOS  curve  for  the  Ti  case.  Since  the  center  of  gravity  of  t 
d  band  and  the  Fermi  level  are  higher  in  energy,  we  should  expect  a  weaker 
interaction  between  the  metal  and  the  CHj.  n  in  Figure  4b  is  pushed  down  less  tha 
in  Figure  4a.  However,  the  main  feature  remains  the  same.  The  Cr  case  is 
intermediate  between  the  Co  and  Ti,  and  we  omit  it. 


CM,  IN  THE  ON -TOP  GEOMETRY 

We  will  analyze  this  specific  case  in  more  detail  than  the  others  so  as  to 
demonstrate  how  the  DOS  and  COOP  curves  may  be  used  to  understand  the  bonding.  Th 
other  geometries  can  be  then  analyzed  in  less  detail.  We  concentrate  on  the  Co 
surface.  The  DOS  curve  of  geometry  7  (Figure  4),  was  already  discussed  above.  We 
have  concluded  that  among  *  •  methyl  orbitals  only  the  n  and  **-CH3  orbitals 
-ontribute  to  metal-carbon  oonding  (i.e.,  interactions^,  ^  and@in  10).  We 
proceed  now  to  examine  these  interactions  in  more  detail. 

We  shall  begin  with  interaction  between  the  n-CH^  orbital  and  the  metal  bands 
If  we  were  in  a  discrete  mono-nuclear  complex  then  we  could  simply  say  that  the  a 
nature  of  the  CHj  n  orbital  allows  interaction  with  dz2,ll,  but  prohibits  them 
with,  say(dxz,  12.  Life  is  not  so  simple  in  the  solid.  Each  metal  orbital 


spreads  out  into  a  band.  Local  interactions  are  dominant,  but  symmetry  limitations 
on  interaction  are  not  so  strong.  We  often  have  to  replace  statements  such  as  "does 
(or  does  not)  interact  with  a  given  level"  by  "interacts  more  (or  less)  with  such 
and  such  part  of  a  band." 
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Let  us  illustrate  this  point  qualitatively  with  the  interaction  of  n.  In  the 
metal  there  is  not  one  dz2  and  dx2 .  orbital  but  many.  13  illustrates 
schematically  some  representative  orbitals  in  the  dz2  and  dx2  bands.  The  orbitals 
at  the  bottom  of  a  band  are  metal-metal  bonding,  those  in  the  middle  non-bonding, 
chose  at  the  top  antibonding.  Things  are  more  complicated  in  three  (or  two) 
dimensions,  but  these  one -dimensional  representations  are  indicative  of  what 
transpires . 


djt  d,f 

13 


What  we  can  say  now  is  that  n  interacts  with  the  entire  dz2  band,  but  perhaps 
more  strongly  with  the  bottom  of  the  band  than  with  the  top,  judging  by  the  overlap 
differences  between  14  and  15.  For  interaction  with  the  dxz  band,  the 
overlap  is  strictly  zero  only  at  the  zone  center  and  edge  (the  most  antibonding  and 
bonding  combinations,  respectively).  It  is  never  very  efficient,  but  as  16 
shows,  one  can  have  an  overlap  between  the  middle  of  the  dxz  band  and  n.  Still,  ch 
overlap  in  16,  depending  as  it  does  on  non-nearest  neighbor  interaction,  is  not 
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Interaction  diagram  17  is  highly  schematic.  Let's  see  the  actual 
manifestations  of  the  bonding  in  the  DOS.  We  already  saw  the  n  peak  move  down  in 
energy  in  the  DOS  decomposition  of  Figure  4.  Now  let  us  examine  the  detailed 
contribution  of  n,  s  and  dz2  to  the  DOS.  This  is  shown  in  Figure  4,  along  with  the 
COOP  curve  of  Figure  6. 

Figures  5  and  6  here 

The  peak  in  the  DOS  of  the  composite  system  at  -13  eV  is  mainly  the  methyl  lone 
pair,  n.  It  is  stabilized  by  interaction  with  surface  dz2  and  s.  This  is  attested 
to  by  the  contribution  of  these  orbitals  --  10%  of  the  total  s  states  and  9%  of  the 
dz2  states  lie  in  this  band  (see  projections  in  Figure  5)  and  their  predominant 
metal-carbon  bonding  character  (see  COOP  curve  in  Figure  6). 

The  antibonding  component  of  the  metal-methyl  interaction  is  also  clearly  seen. 
The  metal  dz2  band,  formerly  confined  to  the  region  of  -12  to  -8  eV,  is  now  broader, 
-13  to  +2  eV.  Much  of  the  density  in  it  is  pushed  up  above  the  Fermi  level.  The 
COOP  curve  shows  a  broad  region  of  M-C  antibonding  from  -11  to  4  eV.  Here  are 
disposed,  highly  delocalized,  the  n-dz2  antibonding  combinations.  The  still  higher 
energy  M-C  antibonding  region  arises  from  out-of-phase  mixing  of  metal  s  and  p,  with 
methyl  lone  pair. 

Some  further  insight  into  the  special  interaction  of  surface  dz2  may  be 
obtained  by  looking  at  the  contribution  to  the  DOS  of  the  surface  Co  atom  not 
involved  in  bonding  to  a  methyl  group,  Co2 ,  and  comparing  it  with  the  bonded  Col. 
This  is  done  in  Figure  7.  Note  how  the  Co2  dz2  DOS  remains  compact,  relatively 
undispersed,  while  that  of  Col  dz2  becomes  quite  spread  out.  Bonding  implies 
dispersion. 


>S  after  the  chemisorption  occurs,  c.  n-CH 


DOS 


DOS 


Figure  7:  A  comparison  of  the  dz2  states,  a  corresponds  to  the  metal  below  the  CH 
group  and  b  the  adjacent  metal  atom  (the  one  not  capped  by  a  CH^) . 


I 


18 

\'e  have  now  seen  how  the  qualitative  model  of  localized  interactions  and  a 
perturbation  theory  based  language  for  discussing  these  interactions  are  beautiful v 
supported  by  the  DOS  and  COOP  curves.  The  molecule  is  bound  to  the  surface 
primarily  by  the  n+dz2,s  interaction.  Many  M-C  antibonding  combinations  are  pus'r.ec 
above  the  Fermi  level.  The  n  band  is  nearly  totally  populated,  with  1.6  electrons 
per  methyl  group.  Note  that  this  brings  the  CH^  group  closer  to  CH-j'  ,  but  then  the 
formalism  of  the  methyl  choice  was  just  that,  a  formalism. 

Another  interaction  that  might  have  led  to  further  charge  transfer  to  CH3  is (7) 
between  various  metal  orbitals,  mainly  dxz ,  d^z  and  methyl  x*-CH3.  The  calculations 
show  clear  signs  of  that  interaction,  but  it  is  not  very  effective,  for  the 
separation  in  energy  between  the  interacting  orbitals  is  large. 

»’e  have  discussed  so  far  the  M-C  bonding,  which  is  obviously  the  major  chemical 
event  that  occurs  when  methyl  is  attached  to  a  metal  surface.  However,  there  are 
additional  bonding  interactions  between  the  metal  and  the  organic  fragment  as  well 
as  changes  within  the  metal  slab  and  within  the  organic  fragment.  Let  us  discuss 
now  briefly  these  "secondary"  changes  which  follow  adsorption. 

Looking  first  at  the  hydrogen-metal  interactions,  we  find  that  these  are 
extremely  small  and  slightly  M-H  antibonding.  This  holds  for  both  the  Ml-H  and  (Ml 
is  the  metal  atom  directly  bonded  to  carbon,  etc.)  and  the  M2-H  interactions.  Thus, 
for  Co,  the  corresponding  overlap  populations  are  -0.022  and  -0.002  respectively.  A 
detailed  examination  of  the  COOP  curve  for  M-H  bonding  shows  an  antibonding 
contribution  in  the  band  at  -12.6  eV,  which  is  mainly  the  bonding  n+dz2,s 
combination.  The  small  M-H  antibonding  effect  is  set  by  the  phase  relationship 
defined  in  18. 
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Our  conclusion  that  the  M-H  interactions  are  weak  and  unimportant  contrasts 
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with  those  of  'fuetterties '  group.  They  have  performed  extended  Hiickel  cluster 
type  calcula’  -.s^a  and  have  reached  the  conclusion  that  multicenter  metal- 
hydrogen-car  interactions  play  a  decisive  role  in  the  chemisorption  of 
hydrocarbons  n  clean  metal  surfaces.  They  emphasized  that  the  most  stable 
geometries  are  those  that  achieve  maximum  multicenter  bonding  of  this  type.  In 
light  of  this  discrepancy  we  looked  further  into  this  problem.  We  changed  the  M-i-H 
angle  from  109.5°  to  95°,  so  that  Gol-H  distance  decreases  from  2.67A  in  6  to 
2 . 45A. 

Still  we  find  no  indication  of  metal -hydrogen  bonding.  On  the  contrary,  as  the 
hydrogens  approach  the  metal  surface  the  antibonding  interactions  increase.  Thus 
our  calculations  predict  repulsion,  not  attraction,  between  the  surface  and  the  ,3- 
hydrogens.  It  is  interesting  that  Minot,  Van  Hove  and  Somorjai^®  used  extended 
Huckel  cluster  calculations  similar  to  those  of  Muetterties  and  did  not  mention  the 
presence  of  strong  metal-0-hydrogen  interactions. 

We  now  turn  to  the  metal -metal  bonding  changes  upon  chemisorption  of  the  methvl 

group.  We  find  that  as  the  new  M-C  bond  is  formed  the  metal -metal  bonds  around  the 

binding  site  are  weakened.  Thus,  the  overlap  population  between  Col  and  Co2  drops 

from  0.185  in  the  bare  metal  to  0.170  in  the  covered  surface  (similarly  for  Cr  and 

Ti  the  values  are  0.462  and  0.411,  0.350  and  0.291,  respectively).  To  some  extent 

the  new  M-C  bond  is  formed  at  the  expense  of  weakening  the  bonds  within  the  metal 

lattice.  For  a  more  detailed  analysis  of  this  general  phenomenon  the  reader  is 

2 1 

referred  to  another  contribution  from  our  group. 

The  overall  effect  of  adsorption  on  the  C-H  bonds  is  small.  The  relevant 
overlap  population  is  0.779  in  isolated  CHj ,  0.787  in  isolated  CH^’  and  0.791  (for 
the  Co  surface). 


Further  insight  into  the  bonding  mechanism  is  provided  by  analyzing  the  charge 
distributions.  The  methyl  is  calculated  to  be  strongly  negative  with  a  charge  of 
0.585  (see  Table  2).  This  is  for  Co.  High  negative  values  are  observed  for  the 
other  metal  surfaces.  Thus,  if  the  metal  and  the  CH^  fragments  are  taken  as  neutrj 


/v 
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when  separated,  then  upon  chemisorption  strong  electron  transfer  from  the  metal  - 
the  methyl  takes  place;  the  metal  becomes  positively,  and  the  methyl  negatively 
charged.  This  is  true  for  Ti  and  Cr  surfaces  as  well  as  Co,  thought  there  are 
differences  to  which  we  will  return  later. 

Table  2  here 

Where  does  this  electron  transfer  come  from,  and  how  can  we  reconcile  it  with, 
the  organometallic  view  of  methyl  as  a  donor?  The  electron  transfer  derives 
primarily  from  the  large  (80%)  occupation  of  n  on  adsorption,  and  that  in  turn  is  a 
consequence  of  the  n+dz2,s  bonding  combinations,  the  large  peak  in  the  DOS  at  -12.-1 
eV,  coming  way  below  the  Fermi  level.  To  put  it  another  way,  and  at  the  same  time 
to  answer  the  second  question  posed  above,  on  the  surface  n  is  filled,  and  a  better 
conceptual  starting  point  might  have  been  to  think  of  it  as  an  adsorbed  CH-j' ,  a 
carbanion.  If  one  views  the  charge  transfer  as  beginning  from  CH^"  and  (slab)+  then 
the  process  of  bonding  is  accompanied  by  a  drift  of  0.415  electrons  from  CH-j*  to  the 
slab.  The  direction  of  this  electron  drift  might  make  some  organometallic  chemists 
happier;  we  think  what  is  important  is  that  one  perceives  the  equally  valid  dual 
viewpoints  of  the  bonding  process. 

From  which  levels  of  the  metal's  "sea  of  electrons"  do  these  transferred 
electrons  come?  We  return  here  to  a  neutral  CH-j  reference  (i.e.  the  calculated 
transfer  being  to  methyl) .  Examination  of  the  electron  distribution  chart  (for  Co 
19  shows  that  most  of  the  electrons  (0.47  electrons)  come  from  one  metal  atom-- 
Col,  which  is  directly  bonded  to  the  organic  fragment.  The  adjacent  atoms  in  the 
surface  layer,  Col,  also  participate  by  donating  0.12  electrons.  The  other  atoms 
are  merely  spectators  in  the  bonding  process.  The  inner  B'and  the  second- inner  C 
layers  remain  essentially  unchanged.  The  only  atom  in  the  inner  layers  that  "feels' 
the  adsorbent  is  the  Co  atom  located  exactly  below  Col,  which  gains  0.07  electrons 
Interaction^  in  10,  which  involves  charge  transfer  within  the  metal  slab,  is 


small  in  this  case. 
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Thus  the  charge  transfer  is  a  localized  phenomenon,  occurring  mainly  between 
the  two  atoms  (Col  and  C)  that  form  the  new  bond.  Furthermore,  not  all  the  orbitals 
on  Col  contribute.  Most  of  the  charge  comes  from  the  dz2  orbital  (0.53  electrons), 
the  dX2  and  d^z  orbitals  contribute  only  0.002  electrons  each.  The  small  electron 
donation  from  the  dxz  and  d^z  orbitals  is  a  strong  indication  that  the  p-type 
interactions  with  are  relatively  weak  and  charge  transfer  occurs  mainly  in 

the  a-framework  through  the  M-n  interaction  (see  17)-  Thus,  based  on  the  charge 
criteria,  M-C  bonding  is  due  mainly  to  the  a  M  +  n  interaction.  Similar  conclusion 
can  be  drawn  for  the  T1  and  Cr  surfaces. 

To  summarize:  the  overall  picture  that  emerges  from  our  analysis  is  that  of 
localized  bonding  where  essentially  only  one  metal  atom  participates  in  the 
formation  of  the  new  metal-CHj  bond.  The  orbitals  that  participate  are  also 
limited.  In  the  metal  slab  these  are  mainly  dz2  and  s  on  Ml  and  on  the  methyl  it  is 
primarily  the  methyl  lone  pair,  n.  The  other  metal  atoms  as  well  as  the  other 
orbitals  hardly  feel  the  adsorbent.  This  bonding  description  is  very  similar  to 
that  in  an  isolated  ML^CH-j  complex,  of  which  many  are  known.  There  is  one 
difference  in  the  surface  case.  In  the  molecular  case  the  ML^  acceptor  orbital  that 
interacts  with  the  CHj  n  is  an  empty  dz2-s-pz  hybrid,  and  that  hybrid  become 
populated  upon  interaction.  In  the  case  of  the  surface  the  orbital  that  interacts 
with  CH^  n  is  not  one  orbital,  but  the  entire  dz2  band.  The  metal  actually  loses 
electron  density  through  this  interaction,  as  some  dz2-n  antibonding  combinations 
are  pushed  above  the  Fermi  level. 


Lee  us  examine  now  how  the  bonding  between  the  methyl  and  rhe  surface  chore¬ 
as  the  electron  count  of  the  meta*  is  decreased.  Relevant  overlap  population  d.it 
for  Co,  Cr  and  Ti  as  well  as  the  corresponding  binding  energies  are  presented  in 
Table  3. 

As  one  sweeps  across  the  transition  series  from  left  to  right  the  average 

energy  of  a  d  electron,  the  center  of  gravity  of  the  d  band,  moves  down  This 

variation  may  have  certain  consequence,  for  as  the  center  of  gravity  of  the  band 

shifts  different  adsorbate  levels  may  interact  more  or  less  strongly  A  good 

illustration  of  this  may  be  found  in  another  reaction  we  have  studied,  CO 
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chemisorption  and  dissociation.  As  one  moves  to  the  left  in  the  transition 
series,  the  CO  n*  becomes  increasingly  occupied,  with  consequent  dissociative 
chemisorption . 

According  to  the  COOP  curve  of  Figure  6,  the  d  band  energy  range  above  the  n 
peak  is  mainly  M-C  nonbonding.  Thus  as  the  d-band  is  depopulated  the  M-C  overlap 
population  (which  is  summed  up  to  the  Fermi  level)  is  expected  to  remain  about  rh 
same.  In  the  Ti  case,  however,  the  d  band  is  higher  in  energy  and  should  ir.terac 
less  effectively  with  n  (by  the  energy  criterion  of  eq.  2).  This  is  reflected  in 
the  M-C  overlap  populations  in  Table  3. 


The  metal -metal  overlap  population  increases  on  going  from  Co  to  Cr,  as 
expected  from  the  fact  that  the  top  of  the  d-band  is  metal-metal  antibonding.  Bu 
from  Cr  to  Ti  some  bonding  states  are  emptied,  so  we  expect  the  overlap  populatio 
to  decrease  and  the  M-M  bond  to  be  weakened.  Figure  8  shows  a  typical  COOP  curve 
for  the  M-M  bond  which  is  the  basis  of  this  reasoning. 


The  binding  energies  of  the  methyl  fragment  to  the  metal  surface,  whic1 
also  given  in  Table  3,  are  naturally  of  special  interest.  The  binding  er.  _s 

defined  here  as  the  energy  of  a  covered  unit  cell  minus  the  combined  energies  of 
neutral  methyl  fragment  and  of  a  unit  cell  of  the  bare  surface.  Although  the 
extended  Hiickel  method  is  not  expected  to  reproduce  quantitatively  the  binding 


TABLE  3.  Bonding  Information  on  the  Surface  as  a  Function  of  the  Metal 
Electron  Count. 


a.  Ei  =  [ E (slab)  +  E (CH3 ) ]  -  E  ; 

Ez  =  [E(slab+)  +  E (CH 3 ~) ]  -  EtQ°  , •  Thus  the  positive  signs  mean  CH3 
or  CH3”  is  bound.  E;  measures  the  covalent  contribution  to  the  total 
binding  energy  Ei. 
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energies  we  do  believe  thac  che  relative  binding  energies  have  chemical 
significance.  Unfortunately,  there  is  no  reliable  experimental  data  which  can  be 
used  to  calibrate  our  calculation. 

We  have  traced  down  above  the  origins  of  the  bonding  interactions  between  the 
metal  and  che  methyl.  The  binding  energy  gives  the  total  energetic  consequences 
thac  occur  as  che  methyl  binds  to  the  surface.  These  include  the  changes  that  occur 
in  all  bonds,  specifically  the  M-C,  M-M,  M-H  and  C-H  bonds.  Why  do  the  binding 
energies  (E^)  that  are  defined  as  E^-[E(slab)+E(CH2) ] '^total  decrease  on  g°ing  from 
Ti  to  Co  although  the  M-C  bonds  get  stronger  and  the  M-M  bonds  alternate  in  the  way 
we  described  above  (Table  3)?  This  apparent  paradox  is  reminiscent  of  a  similar 
situation  which  occurs  in  using  frontier  orbital  arguments  in  analyzing  the  strength 
of  simple  two  electron  bonds.  Within  this  simplistic  theoretical  framework  it 
seems  chat  the  strength  of  a  bond  should  increase  as  the  energy  separation  between 
the  interacting  orbitals  decreases,  i.e.,  as  its  covalent  character  increases. 

Highly  ionic  bonds,  where  the  separation  between  the  orbitals  that  form  the  bonds  is 
very  large,  are  therefore  expected  to  be  extremely  weak,  contrary  to  reality.  As 
we  use  here  essentially  frontier  orbital  logic  (eq.  2)  in  analyzing  the  M-CH^  bond, 
which  on  the  basis  of  the  charge  distribution  is  substantially  ionic,  it  is  not 
surprising  that  similar  problems  arise.  Picture  20  describes  schematically  the 


energetic  consequences  of  the  ionic  character  of  the  M-C  bond.  The  M-C  bond 
strength  arises  mainly  from  the  fact  that  the  approaching  CH^  radical  brings  along 
a  "hole"  at  an  energy  of  -11.75  eV,  3.26  eV  below  the  Fermi  level  of  cobalt.  The  M- 
C  bond  is  formed  by  an  electron  transfer  from  the  metal  into  this  "hole"  and  a 
binding  energy  of  3.26  eV  is  gained  in  this  process.  Thus,  the  "ionic  character"  of 
the  M-C  bond  can  be  described  as  responsible  for  87%  of  the  total  binding  energy 
(3.73  eV)  of  the  methyl  radical  to  the  Co  surface.  The  difference  of  0.47  eV  can  be 
described  as  the  "covalent  contribution"  to  M-C  bonding  (E2  value  in  Table  3,  see 
also  discussion  below) . 
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As  we  move  from  Co  to  Ti  the  Fermi  level  changes  by  -2  eV  to  higher  energy. 

The  energy  gain  on  filling  the  "hole"  is  now  greater  (see  20) ,  and  so  is  the 
binding  energy  and  charge  transfer  to  the  methyl  (-0.59  and  -0.79  for  Cr  and  Ti 
respectively).  The  smaller  "ionic  contribution"  to  the  binding  energy  for  Co 
compared  to  Ti  is  partially  compensated  by  the  increased  "covalent  character"  of  the 
Co-C  bond  which  is  reflected  in  the  increased  M-C  overlap  population  (Table  3).  The 
binding  energy  of  methyl  to  Cr  is  therefore  lower  than  that  of  Co-CH3  by  only  1.7 
eV,  not  2  eV  as  suggested  in  20. 

Once  again  it  is  useful  to  change  perspective  and  think  of  the  binding  energy 
of  a  CH-j"  anion  to  a  positively  charged  slab.  The  "hole"  is  filled  and  ambiguities 
of  charge  transfer  are  avoided.  Now  the  binding  energy  E2 ,  defined  as 
E2-iE(slab+)+E(CH3') ] -Etotal,  is  lower  for  Ti  (0.16  eV)  than  for  Co  (0.47eV), 
paralleling  the  M-C  overlap  populations.  Cr  is  intermediate  between  Ti  and  Co;  due 
to  its  nonhexagonal  lattice  it  might  deviate  from  average  behaviour. 

A  final  point  to  be  made  about  the  on-top  site  is  that  the  barrier  for  the 
rotation  of  the  methyl  group  around  its  local  C3  axis  of  symmetry  (see  21)  is 
nearly  zero.  This  is  hardly  surprising,  for  this  is  a  sixfold  barrier,  and  we  know 
such  to  be  very  small  in  molecular  cases. 


Figure  9  displays  the  total  DOS  (dashed  line)  of  the  Co(0001)  system  upon 
chemisorption  of  CH^  in  a  doubly  bridging  geometry.  The  darkened  areas  show  the 
contributions  of  the  organic  fragment  and  the  arrows  indicate  the  location  of  thes< 
states  in  the  isolated  fragment.  Figure  10  shows  the  COOP  curves  for  the  metal- 
carbon  bonds.  Comparison  of  Figures  9  and  5  and  of  10  and  6  reveals  immediately  c 
similarity  in  the  bonding  mechanisms  in  the  on- top  and  the  bridged  geometries.  Ir. 
both  cases  the  interaction  which  is  primarily  responsible  for  the  metal-carbon 
bonding  is  that  with  the  n  orbital.  Nevertheless,  there  are  important  differences 
in  the  details  of  the  bonding  mechanisms  in  the  two  geometries  and  we  concentrate 
our  attention  on  these. 

Figure  9  and  10  here 


The  major  difference  is  obvious.  In  the  "on- top"  geometry  only  one  surface 
atom  contributes  to  M-CH2  bonding,  while  in  the  bridging  case  two  metal  atoms  are 
clearly  involved.  More  interesting  is  the  contribution  of  the  various  orbitals  to 
M-C  bonding.  In  contrast  to  7,  the  carbon  atom  in  8  is  located  off  the  dxz 
nodal  plane  and  effective  Interaction  with  the  n  orbital  can  now  occur,  see  22 
(contrast  with  12) . 
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As  seen  in  22  these  interactions  take  place  with  dxz  levels  at  or  near  the 
edge  of  the  Brillouin  zone,  i.e.  the  top  of  the  band.  Projection  of  the  dxz  leve, 


(Figure  11)  indeed  shows  considerable  broadening  of  this  band  compared  with  6 


Figure  9:  Total  DOS  (dashed  line)  and  the  CH3  contribution  (darkened  area)  of  the 
bridging  CH3  +  Co(0001)  system.  The  arrows  indicate  where  the  free  CH-j  MO  levels 
were  before  the  adsorption. 


les  are  the  integr 


Figure  5).  In  the  bridging  geometry  the  dxz  band  contributes  12%  (6%  for  each 

the  Co  atoms)  of  the  total  number  of  states  in  the  n  band  at  -13ev.  In  7  this 
contribution  was  close  to  zero.  Another  consequence  of  bridging  is  that  n  interac 
effectively  only  with  the  bottom  of  the  dz2  band  (23),  and  not  with  the  top 

2-0  . 


We  thus  expect  that  the  contribution  to  M-C  bonding  from  the  dz2  band  will  be 
smaller  in  the  bridging  mode  than  in  the  on-top.  This  is  indeed  observed.  The 
contribution  of  the  dz2  band  to  M-C  bonding  (i.e.  its  portion  of  the  n  peak  in  the 
DOS)  drops  from  9%  in  7  to  less  than  2%  in  8  (compare  Figures  11  and  6a).  The 
participation  of  the  s-band  also  diminishes  from  16%  in  7  to  12%  in  8. 

The  contribution  of  the  w-type  interaction  to  M-C  bonding  is  higher  in  the 
bridging  than  the  on- top  site.  This  is  indicated  in  Figure  10  by  the  increased 
contribution  of  the  x-CH^  orbital  at  -15.3  eV  to  M-C  bonding  (compare  with  Figure 
7).  A  major  difference  between  8  and  7  is  that  in  the  former  the  degeneracy  of 
the  methyl  n- interactions  is  destroyed.  These  are  two  distinct  components,  one 
pointing  along  the  M-M  axis  (i.e.,  x  in  25)  and  the  other  in  a  direction 
perpendicular  to  it  (y  in  26).  The  *x  orbital  interacts  more  strongly  with  the 
metal  surface  orbitals  (see  25)  than  the  ity  orbital  (see  26).  This  is  due 
primarily  to  the  poorer  orbital  overlap  in  26,  but  also  to  the  energy  of  the 
piece  of  the  dxz  or  dyz  band  with  which  overlap  is  effective.  The  representative 
dxz  orbital  is  near  the  top  of  its  band  (see  the  nodal  relationship  in  25)  and 
so  closer  in  energy  to  its  w*-CH.j  partner.  The  representative  d^z  orbital  is  low 


its  band,  26. 
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This  effect  is  beautifully  exhibited  by  comparing  Figures  11c  and  lid  which 


ft  ft 

'.ow  the  projected  DOS  curves  of  the  methyls  n  and  n  orbitals  respectively. 

x  y 


Vhile  the  n  orbital  is  dispersed  into  a  band  of  significant  width  indicating 


strong  interaction  with  the  metal,  the  »  ^  band  remains  relatively  narrow. 


The  two  bonded  metal  atoms  contribute  similarly,  as  expected,  to  the  bonding 
the  CHj  fragment.  This  is  apparent  from  the  total  M-C  overlap  populations^  shown  : 
2/that  indicate  a  slightly  stronger  Col-C  bond.  Col  is  "staggered"  with 
respect  to  the  methyl  hydrogens,  Co2  "eclipsed". 
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Closer  examination  of  the  COOP  curves  in  Figure  10  reveals  an  intriguing 
difference  between  Col  and  Co2.  The  COOP  peak  corresponding  to  the  M+n  interact! 
is  almost  twice  as  large  for  Co2  than  for  Col.  In  the  d-band  region  the 
contribution  of  Co2  to  M-C  antibonding  is  also  larger  than  that  of  Col.  For  coba 
where  the  d-band  is  almost  filled,  the  stronger  contributions  of  Co2  to  M-C  bondi 
and  antibonding  nearly  compensate  and  the  Col-C  and  Co2-C  overlap  populations  are 
similar,  the  former  being  somewhat  larger.  As  the  d-band  is  depopulated  the  M2-C 
bond  should  gradually  strengthen.  The  overlap  populations  calculated  are  in  fact 
reversed  for  M  -  Ti,  compared  to  Co  (see  27). 

The  stronger  M2-C  Interaction  In  particular  in  the  n  band^ must  result  from  t 
|  different  arrangement  of  the  methyl  hydrogens  with  respect  to  the  two  relevant  me 

atoms.  There  is  a  correlated  differential  in  the  direct  M-H  and  C-H  overlap 

I 


populations  (see  28  and  29) .  But  the  effect  is  small  and  we  are  not  certain 


if  it  merits  a  detailed  discussion. 
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The  total  interaction  between  the  methyl  fragment  and  the  metal  surface  seems 
stronger  in  the  bridged  geometry  than  in  the  on- top  geometry.  Thus,  in  8  the  M+r. 
band  that  is  the  major  contributor  to  M-C  bonding  is  pushed  to  a  lower  energy  by  i  . 
eV  (for  Co)  relative  to  its  position  in  7  -  indicating  a  stronger  bonding 
interaction.  Also,  the  M-C  overlap  population  data  indicates  overall  stronger  M- 3 
bonding  in  8. 

Although  the  bonding  between  the  metal  atoms  and  the  methyl  is  stronger  in  the 
bridged  geometry,  the  binding  energy  of  the  organic  fragment  to  the  surface  is 

lower  in  8  compared  to  7  see  Table  4).  Thus,  for  Co,  Cr  and  Ti  the  binding 

I 

energies  of  methyl  in  the  bridging  geometry  are  by  1.1,  0.9  and  0.5  eV  respectively 
lower  than  in  the  on- top  site.  This  apparent  paradox  of  a  stronger  M-C  bond  that 
j  results  in  a  weaker  binding  of  the  CH^  fragment  to  the  surface  becomes 

<  comprehensible  upon  examining  the  M-M  overlap  populations  in  30. 
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30  shows  that  the  M1-M2  bond  is  considerably  weakened  on  going  from  7  to 
8.  This  bond  weakening  is  a  direct  consequence  of  the  stronger  metal-methyl 
interactions  that  push  many  of  the  metal  bands  to  higher  energies.  As  the  top  of 
the  d-band  is  metal-metal  antibonding  (see  Figure  8,  the  COOP  for  7  is  similar). 
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these  bonds  are  weakened  as  che  metal - fragment  interactions  become  stronger.  7h<- 
dominant  n-d  interaction  pushes  up  above  the  Fermi  level  those  metal  states  which 
contribute  to  either  M-M  bonding  or  M-M  antibonding  in  the  on- top  geometry.  But 
only  M-M  bonding  states  are  lifted  above  the  Fermi  level  by  this  interaction 
at  the  bridging  position,  and  the  M-M  bond  is  weakened  more.  Thus,  the  stronger  M- 
bonding  in  8  is  gained  in  the  expense  of  weakening  the  bonds  within  the  metal 
slab,  therefore  raising  its  energy.  In  the  Co  case  (and  to  a  smaller  degree  also 
for  Cr  an  Ti)  the  M-M  bonding  energy  that  is  lost  is  greater  than  the  M-C  bonding 
energy  that  is  gained  upon  bridging  and  the  binding  energy  of  methyl  is  lower  in  8 
than  7.  In  addition,  interaction  4  in  10  plays  a  role:  as  the  M-M  bonding 
states  are  lifted  up  in  energy,  the  electrons  flow  from  these  states  into  the  top  o 
the  d-band,  which  is  M-M  antibonding.  Consequently,  when  the  electron  count  on  the 
metal  decreases  and  the  M-M  antibonding  states  are  depopulated  the  difference  in 
binding  energy  between  the  on- top  and  the  bridging  geometry  is  reduced.  In  our 
calculation  the  difference  is  1.1  eV  for  Co  and  0.9  eV  for  Cr,  the  on-top  position 
being  more  stable. 

Only  a  few  molecular  complexes  are  known  that  have  a  methyl  group  bridging  tw 
bonded  metal  atoms^  (e.g.  31^c  and  32^).  In  most  such  complexes  the 
methyl  bridges  two  non-bonded  metal  atoms,  as  in  33.  Thus,  it  seems  that  in 
analogy  to  our  computations  for  the  covered  surface,  in  the  molecular  complexes  the 
bridging  mode  is  accompanied  by  M-M  bond  weakening  or  by  M-M  bond  cleavage. 


31  32  33 


Finally,  we  comment  on  the  charge  transfer  that  occurs  upon  bridging.  This  is 
shown  in  34  and  35  for  the  cobalt  surface. 


34  35 


As  we  can  see  the  total  transfer  of  charge  from  the  metal  to  the  methyl 
fragment  is  nearly  the  same  in  7  and  8  (0.585  and  0.531  electrons  respectively). 
However,  there  appear  to  be  dramatic  changes  in  the  charge  distributions  within  th 
metal  slab  on  going  from  7  to  8 .  There  is  a  substantial  flow  of  charge  from  the 
surface  layer  into  the  inner  layers  upon  methyl  bridging.  The  same  is  observed  fo 
Cr  and  Ti  surfaces.  Thus  for  cobalt  the  calculations  show  that  on  going  from  7  to 
8  the  surface  A'  layer  loses  0.52  electrons  while  the  B'  and  C’  layers  gain  0.30 
and  0.36  electrons  respectively.  Thus,  the  Co  surface  layer  which  in  the  bare  met 
is  negatively  charged  becomes  nearly  neutral  for  on  top  bonding,  and  turns 
positively  charged  in  the  bridging  case.  This  substantial  charge  flow  away  from  t 
surface  is  easily  understood  in  terms  of  the  bonding  picture  that  emerged  from  our 
calculations.  Bridging  involves  stronger  interactions  with  the  metal  d-bands.  Th 
surface  bands  that  are  involved  in  these  interactions  are  pushed  up  to  higher 
energies  and  are  therefore  emptied.  Charge  flows  to  lower  metal  bands  which  are 
concentrated  in  the  inner  slab  layers .  These  charge  flows  provide  a  clear 
determination  of  the  importance  of  interaction (Z) in  10. 

As  in  the  on-top  geometry,  so  also  in  the  bridging  site  the  methyl  is 
essentially  freely  rotating;  there  is  no  calculated  rotational  barrier. 


Once  we  have  analyzed  in  detail  the  on- top  and  the  bridged  geometries  and  hav, 
understood  the  bonding  mechanism,  then  it  is  easy  to  predict  what  happens  in  the 
capping  geometry  9.  The  trends  that  have  been  observed  on  going  from  7  to  8 
will  hold  also  for  moving  from  8  to  9 .  In  9  the  interactions  between  the 
methyl  and  the  metal  slab  are  the  strongest,  and  the  total  M-C  bonding  increases,  a? 
expected,  on  going  from  8  to  9 .  36  and  37  show  the  overlap  populations 

for  the  M-C  and  M-M  bonds  in  both  geometries.  Concurrently  on  going  from  8  to  9 
the  M-M  bonds  are  weakened  as  reflected  in  the  decreasing  M-M  total  overlap 
populations  (see  36  and  37).  At  high  electron  counts  this  results  in  a 
smaller  binding  energy  for  the  capping  geometry  relative  to  the  on-top  mode  (e.g., 
by  0.3  eV  for  Cr) .  For  lower  electron  counts,  many  of  the  M-M  bonding  states,  whici 
are  pushed  up  above  the  Fermi  level  by  the  CH^-metal  interaction,  were  originally 
empty.  The  gain  in  M-C  bonding  thus  compensates  more  the  loss  in  the  M-M  bonding  i: 
the  lower  electron  count  case.  On  the  Cr  surface  our  calculation  shows  no  energetic 
difference  between  the  bridging  and  capping  sites  and  on  Ti  the  capping  geometry  is 
more  stable  by  0.3eV.  In  all  the  cases  under  consideration,  however,  the  on- top 
mode  is  the  most  stable  one  for  CH-j . 
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The  charge  distributions  also  behave  as  expected.  There  is  substantial 
electron  flow  from  the  surface  layer  into  the  inner  metal  layer.  The  charge  on  the 
methyl  fragment  decreases  as  its  coordination  number  increases  on  the  Co  surface. 
Thus  the  CH^  charge  is  -0.58  in  7,  -0.53  in  8  and  -0.49  in  9.  This  trend  can 
also  be  rationalized  within  our  bonding  model  using  the  following  reasoning:  If 

there  is  no  interation  between  the  metal  and  the  CH2  radical  and  the  electrons  are 
assigned  to  the  lowest  available  levels,  then  the  methyl's  n  orbital  which  is  lover 


in  energy  chan  che  mecal's  d-band  will  be  doubly  occupied  and  che  charge 
distribution  is:  C^’^metal  slab]  +  .  When  the  M-CH^  interactions  are  "turned  on", 
part  of  the  methyl  charge  is  transferred  back  to  the  metal.  As  the  number  and 
efficiency  of  the  interaction  mechanisms  increases,  as  is  the  case  along  the 
progression  7  -->8  -->9,  more  charge  is  transferred  back  to  the  metal  and 
the  organic  fragment  becomes  less  negative. 

The  capping  or  triply  bridging  geometry  is  the  first  where  we  find  a 

significant  barrier  of  6.3,  21.2  and  21.7  kcal/mol  (for  Co,  Cr,  Ti)  for  the  rotation 

of  the  methyl  group.  The  most  stable  arrangement  is  such  that  the  hvdrogens  point 

towards  three  metal  atoms.  A  side  view  is  shown  in  38.  The  less  stable 

arrangement  in  which  the  hydrogens  point  towards  the  center  of  three  M-M  bonds  is 

shown  in  39.  The  M-H  overlap  populations  which  are  given  in  38  and  39 

suggest  that  weak  M-H  attractive  interactions  favor  38  over  39,  where  these 

interactions  are  repulsive.  38  is  the  only  case  where  we  find  some  indication 

for  the  attractive  M-H  forces  that  Muetterties  and  co-workers  found  to  be 
22 

significant.  The  larger  rotational  barrier  at  lower  electron  count  is 

due  to  the  following  fact:  At  the  eclipsed  geometry  38  the  hydrogens  of  the  CH3 

group  interact  with  both  M-M  bonding  and  antibonding  states,  whereas  in  the 

staggered  39  they  interact  only  with  the  M-M  bonding  states.  Thus  38  is 

more  stable.  At  lower  band  filling  it  is  more  so  because  the  unfilled  M-M 

antibonding  states  act  as  acceptors.  However,  we  are  not  sure  whether  the  extended 

Huckel  method  has  overestimated  the  barrier.  Better  calculations  are  needed. 
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The  computational  results  for  all  three  sites  of  methyl  binding  are  summarized 


in  Table  4. 
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To  conclude  chis  section  we  summarize  the  main  features  of  the  bonding  be  tv 
3.  methyl  fragment  and  the  metal  slab.  The  major  interactions  occur  between  the 
metal's  surface  bands  and  the  CHj  orbital.  These  interactions  are  a  type,  invol 
surface  d,2  and  s  orbitals.  The  important  difference  from  a  molecular  complex  i 
that  charge  transfer  in  the  a  system  occurs  from  d22  (and  s)  to  CHj.  But  this  i 
consequence  of  the  reference  state  as  neutral  CH-j ,  and  if  CH-j'  is  chosen  instead 
the  methyl  group  behaves  in  the  way  we  would  have  expected.  As  one  moves  along 
geometry  series:  on-top  -->  bridged  -->  triply  bridging  some  additional  M-C  bond 
is  gained  but  this  is  generally  achieved  at  the  expense  of  considerable  weakenin 
the  metal-metal  bonds  in  the  surface  layer.  As  a  result  the  on-top  geometry  is 
most  stable  geometry  for  all  three  metal  surfaces. 


ADSORBED  METHYLENE 

k’e  will  follow  the  guidelines  of  our  analysis  for  the  methyl  case.  As  in  t 
methyl  system  three  adsorption  sites  are  studied:  the  1-fold  on-top  geometry 
jO,  the  2-fold  bridging  position  41  and  the  3-fold  bridging  or  capping 
position  42.  For  each  of  these  geometries  two  conformers  with  respect  to 
rotation  around  the  C2V  axis  of  the  methylene  group  were  examined. 
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The  major  difference  in  the  bonding  of  methyl  and  of  methylene  to  the  metal 
surface  lies  on  the  different  set  of  orbitals  that  these  two  fragments  bring  to 
bonding  process.  The  relevant  orbitals  of  methylene  are  displayed  in  43  in 
increasing  order  of  energy  (in  parenthesis,  eV)  from  left  to  right.  In  the  free 


,  ir.glet  methylene  only  and  a  are  occupied.  Note  that  for  consistency  vith 

-ethyl  case  we  should  have  called  the  lone  pair  orbital  n,  but  it  is  customary  t: 


iesignate  it  as  a  in  carbene  chemistry. 
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A  schematic  energy  diagram  of  the  methylene  orbitals  and  the  metal  slab  leve 
Ls  displayed  in  44.  44  also  shows  the  position  of  the  relevant  orbitals  of 
:he  methyl  radical. 
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Let  us  examine  first  the  most  important  features  of  the  surface  -  carbene 
oonding,  those  which  hold  for  all  the  binding  sites.  Diagram  44  reveals  at  a 
glance  the  important  differences  in  the  bonding  of  methyl  and  methylene  to  the 
surface.  The  methyl  fragment  carries  a  singly  occupied  empty  orbital  at  -11.75eV 
while  the  methylene  introduces  an  entirely  vacant  orbital  at  a  similar  energy, 
have  seen  above  that  most  of  the  binding  energy  of  the  organic  fragment  to  the  me 
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surface  can  be  associated  with  the  "ionic  character"  of  che  M-C  bond,  in  other 
words  with  the  energy  gained  upon  filling  the  "hole"  that  the  organic  fragment  is 
carrying.  Using  this  naive  picture  we  expect  that  methylene  which  carries  2  emrtv 
sites  should  be  bound  much  more  strongly  to  the  surface  than  methyl,  which  proviso, 
only  one  vacancy.  Note  also  that  the  empty  orbital  that  methylene  brings  along  is 
of  p  or  k -  type  while  that  of  methyl  is  of  a- type.  The  bonding  of  methylene  to  the 
metal  surface  is  therefore  expected  to  have  a  much  more  pronounced  rr-character  th .it 
for  methyl.  We  will  see  below  that  these  qualitative  expectations  are  indeed 
fulfilled. 

Another  important  difference  between  methyl  and  methylene  is  that  in  CH0  the 
r*-orbital  is  doubly  occupied  while  in  CH-j  it  is  singly  occupied.  This  orbital 
interacts  mainly  with  filled  metal  levels  (see  44).  When  two  filled  orbitals 
interact  this  results  in  a  strong  destabilization  due  to  the  occupation  of 
antibonding  orbitals.  This  would  have  indeed  been  the  case  for  the  interaction 
of  the  cr-orbital  of  CH2  with  a  filled  orbital  of  another  molecular  fragment. 

However,  the  situation  in  the  metal  slab  is  entirely  different.  The  electrons  whir:, 
should  enter  the  antibonding  states  can  instead  occupy  empty  metal  levels  which  are 
lying  near  the  Fermi  level.  These  states,  which  we  have  shown  above  to  be  primarily 
bulk- like,  serve  as  an  empty  reservoir  for  electrons  that  are  pushed  up  by 
interactions  between  filled  levels.  These  arguments  are  displayed  schematically  in 
-o . 
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The  strong  desta.  c  that  Is  associated  with  4-electron  interactions  i 

.r.olecules ,  and  that  ccr.  r.uch  of  their  chemistry,  has  to  a  large  extent 

diminished  in  metal-adsc._  ..t  systems.  This  is,  we  believe,  one  of  the  importar. 

cnaracter tstics  of  metals  and  metal  surfaces,  a  factor  that  makes  them  behave  so 

differently  from  the  "analogous"  isolated  metal  complexes. 

We  can  proceed  now  to  discuss  in  detail  the  different  binding  sites  of 

methylene.  Table  5  and  6  collects  the  computational  data  regarding  the  bonding 

characteristics  of  methylene  in  the  different  adsorption  sites.  Some  of  this  dat 

will  be  discussed  in  some  detail  in  the  text.  We  also  refer  the  interested  reade 

to  a  recent  contribution  from  our  group  in  which  we  studied  the  adsorption  of 

2 1 

vinylidene,  HjC-C:  on  metal  surfaces.  In  many  aspects  this  system  is  similar  t 
ours  and  we  will  refer  to  it  when  appropriate. 


We  have  tied  two  conformers  shown  in  46  and  47  differing  by  a  30° 
rotation  arour.  .  the  metal-carbon  bond.  We  find  that  the  energy  difference  betvee 
these  conformers  is  essentially  zero.  We  will  comment  on  the  reasons  for  this  ti 
barrier  later.  We  concentrate  therefore  on  one  geometry,  46,  using  the 
indicated  coordinate  system  (but  the  analysis  applies  to  47  as  well. 
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Figure  12a  and  12b  shows  the  total  DOS  curve  (dashed  line)  of  the  Co  system 
before  and  after  the  chemisorption  of  CH2  in  this  geometry.  The  dark  area  in  Fig. 
12b  gives  the  states  contributed  by  the  methylene  fragment  and  the  straight  lines 
Figure  12a  indicate  the  location  of  these  states  in  the  isolated  organic  fragment 
Figure  12b  can  be  compared  with  the  analogous  Figure  4  for  the  1-fold  chemisorbed 
methyl.  The  two  orbitals  of  methylene  that  are  most  affected  by  the  interaction 
with  the  metal  are  a  and  p.  These  are  pushed  to  lower  energies  by  1 . 2eV  and  O.-e'. 
respectively.  The  peak  of  the  cr-band  is  located  at  a  slightly  lower  energy  in 
2 6  than  in  7.  This  is  due  primarily  to  the  fact  that  er-Ct^  is  lower  than 
n-CH^  to  start  with  (-12.1  vs.  -11.8eV). 


Figure  12...here 

The  fact  that  the  a  orbital  is  affected  more  strongly  than  the  p  orbital  by 
interaction  with  the  metal  indicates  that  the  M-C  bond  is  mainly  of  a  character. 
This  interpretation  is  corroborated  by  Figure  13  which  displays  the  COOP  curve  for 
the  Col - C  bond.  The  strongest  M-C  bonding  peak  is  at  -13.0eV  corresponding  to  the 
a  orbital  of  methylene.  The  p  orbital  also  makes  an  important  contribution  to  M-C 
bonding,  as  indicated  by  the  peak  at  -11.2eV.  But  it  does  so  to  a  smaller  degree 
than  the  a  band. 

Elams  11  here 

Comparison  of  Figure  13  with  the  analogous  COOP  curve  for  methyl  (Figure  6)  is 
interesting.  There  are  essential  similarities,  and  small  differences  -  the  absence 

of  a  p  bonding  contribution  in  Figure  6,  minor  and  understandable  differences  in  the- 

9 

contribution  of  the  *  orbitals.  Note  also  that  around  the  Fermi  levels  of  Co(d  ) 
indicated  by  the  arrow  in  Figure  13,  the  M-C  COOP  curve  is  antibonding  and  it  is 
more  so  for  methylene  than  in  the  methyl  system.  Thus  as  the  top  of  the  d-band  is 
depopulated  e  M-C  bond  is  strengthened.  The  M-C  overlap  populations  are  in 
agreement  .,.441  for  Co  and  0.520  for  Cr  (Table  5).  As  in  the  methyl  case,  for 
r.ethylene  the  M-M  COOP  curve  (which  is  also  shown  in  Figure  13)  is  also  antibondirg 
3round  the  Fermi  level  of  Co.  Emptying  the  top  of  the  metal  band  results  therefor-; 


bonded  to  the  CH,  group  (magnified). 


Figure  13:  COOP  curves  of  the  M-M  (solid  line)  and  the  C-K  (dotted  line)  bonds  fo 
the  on- top  CH2  +  Co(0001)  systea. 


ir.  a  strengthening  of  the  M-M  bond,  up  to  a  certain  point  (-  d5),  past  which  it  :  . 
weakened.  The  calculated  overlap  populations  for  the  on-top  methyl  and  methvler.e 
geometries  are  shown  in  48  (see  also  Table  5).  Note  that  also  for  methylene  the 
M-H  interactions  are  weakly  antibonding. 
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The  calculated  total  overlap  populations  in  48  indicate  that  the  M-C  bond 
is  stronger  for  methylene  than  for  methyl.  This  comes  from  from  the  increased 
contribution  of  ir-bonding.  The  difference  in  the  overlap  populations  between  the 
two  cases  is,  however,  too  small  to  indicat^the  presence  of  a  fully  developed  M-C 
double  bond  character.  Qualitatively,  the  M-CH2  bond  is  best  described  as  a  parria 
double  bond. 

The  relatively  modest  increase  in  the  M-C  overlap  population,  and  thus  in  the 
M-C  bond  strength,  upon  substitution  of  methyl  by  methylene  does  not  indicate,  as  . 
have  already  emphasized  earlier,  that  the  binding  energies  of  these  fragments  to  th 
metal  surface  are  similar.  This  is  because  the  overlap  populations  reflect  only 
what  we  have  termed  the  "covalent"  fraction  of  the  bonding  (E2)  while  most  of  the 
binding  energy  comes  from  an  "ionic"  transfer  of  electrons  from  the  metal  to  the 
organic  fragment.  For  the  adsorption  of  CH2  on  cobalt  we  calculate  a  binding  enert 
of  6.32eV  compared  with  only  3.73eV  for  methyl.  These  values  agree  well  with  our 
interpretation.  The  "ionic"  binding  energy  is  given  by  the  energy  difference 
between  the  top  of  the  bare  cobalt  band  and  the  a-orbital  "hole”,  multiplied  by  2 
(for  2  electrons),  i.e.,  2x(ll .40-8 .48)  -  5.84eV.  The  remaining  fraction  of  the 
binding  energy:  6.32-5.84  -  0.48eV,  can  be  associated  with  the  "covalent"  bonding. 
For  methyl,  only  one  electron  contributes  to  the  "ionic"  stabilization;  i.e., 


.L.^5-3.-%3  -  3.2/eV,  and  the  total  binding  energy  is  only  3.73eV  (leaving  0.46eV  : 

' covalent"  stabilization),  roughly  half  than  that  of  methylene. 

On  going  from  Co  to  Cr  to  Ti,  the  Fermi  level  becomes  higher  and  the  "ionic'1 
energy  gain  is  larger,  increasing  the  total  surface -adsorbent  bonding  interact i or.  - 
see  values  in  Table  5).  The  "covalent"  bonding  remains,  as  expected,  roughlv 
the  same  as  the  metal  is  altered  (see  E2  values  in  Table  5). 

Now  that  we  have  discussed  the  general  features  of  the  bonding  of  CH0  to  the 
cecal  surface  let  us  analyze  the  bonding  of  CH2  to  the  metal  surface  in  some  more 
detail.  We  have  shown  using  Figures  12  and  13  that  M-C  bonding  is  associated  mai: 
with  the  a  and  p  orbitals  of  CH2 .  Which  of  the  metal  levels  contribute  more  to  c: 
M-C  bond? 

Interactions  with  the  n  orbital  of  CH2  involve  mainly  the  metal's  s  and  d„2 
bands  as  shown  in  49  and  50  respectively.  Interaction  49  and  50  are 
analogous  to  the  a  interactions  with  the  methyl  fragment  shown  earlier  in  13. 
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As  for  adsorbed  CH3,  so  also  for  CH2  essentially  the  entire  s  and  d22  bands 

involved  in  the  interaction.  This  is  the  reason  that  we  don't  show  the  orbital 

phases  at  the  neighboring  metal  atom.  Projection  of  the  dz2  states  of  Col  shown 

Figure  12c  reveals  that  nearly  8%  of  these  states  contribute  to  the  metal-<7  band. 

For  the  s  bands  (Figure  12d)  the  contribution  is  slightly  larger  (12%).  The  stro: 

interaction  of  the  dz2  band  with  the  organic  fragment  is  also  evident  from  its  la 

dispersion;  nearly  40%  of  the  dz2-states  are  above  the  Fermi  level  of  cobalt. 

The  p  orbital  of  CHj  can  interact  with  several  metal  bands.  The  major 

interaction  is  with  the  d  levels  as  shown  in  51.  According  to  Figure  12e,  17% 

yz 

of  the  d.,z  states  are  involved  in  w- interactions  with  the  p  orbital.  The 


neighboring  metal  atoms  also  contribute,  as  emphasized  in  51  by  the  arrows.  The 


40 

4..,  levels,  being  perpendicular  to  the  p  orbital,  do  not  interact  with  the  organ:, 

•ft 

fragment  (except  for  a  minor  interaction  with  the  it  -CH2  orbital/ which  is  of  prop-.-, 
svmmetrv  but  lies  18eV  above  the  p  orbital).  Some  levels  around  the  center  of  the 
4,2  and  s  bands  can  also  contribute  to  it  bonding,  as  shown  for  dz2  in  52.  These 
interactions,  however,  do  not  in  fact  contribute  much  to  the  total  it  bonding, 
because  of  poor  orbital  overlap.  Only  1-2%  of  the  dz2  levels  contribute  to  the  M-t 
3  and . 


We  can  now  understand  why  the  barrier  to  the  rotation  of  the  CH2  group  is  so 
small.  The  a  type  interactions  are  invariant  to  rotation  of  the  CH2  group  and  are 
the  same  in  46  and  47.  The  important  interactions  for  the  rotation  are  the 
t  type.  As  the  CH2  fragment  is  rotated  towards  geometry  47  the  contribution  of 
the  d,/z  bands  to  it  bonding  decreases  due  to  smaller  orbital  overlaps.  However,  at 
the  same  time  the  contribution  of  the  dxz  bands  (which  are  zero  in  46)  increases 
ar.d  reaches  its  maximum  in  geometry  47.  A  nearly  constant  contribution  of  these 
levels  to  the  M-C  bonding  is  maintained  as  the  CH2  fragment  rotates.  The  rotation 
barrier  is  six-fold  and  tiny.  The  different  effectiveness  of  the  a  and  it 
Interactions  is  manifested  also  in  the  appropriate  metal  DOS. 

In  55,  we  look  at  the  charge  distribution  in  the  adsorbed  methylene.  We 
find  two  comparisons  to  be  useful,  one  with  the  bare  metal  54  and  the  other  with 
the  analogous  1-fold  on- top  methyl  system  53.  The  data  in  53-55  is  for 
the  Co  surface  and  data  for  other  electron  counts  can  be  found  in  Table  5.  As  the 
CH-,  fragment  provides  2  vacancies  for  the  metal  electrons  while  the  methyl  provide 


only  one,  the  charge  transfer  to  the  organic  fragment  is  much  higher  in  the 


Charges 


methylene  system.  The  CH2  charge  is  -1.54  compared  with  -0.58  for  CH2 .  Using  the 

adsorbed  methyl  system  as  the  reference  point  we  calculate  that  abstraction  of  a 

hydrogen  atom  causes  0.90  electrons  to  flow  from  the  metal  to  the  organic  fragment. 

The  charge  is  donated  roughly  equally  from  all  three  metal  layers.  Comparison  with 

the  bare  surface  shows  that  as  in  the  methyl  system,  the  charge  that  is  donated  to 

the  methylene  comes  mainly  from  the  Co  atom  that  is  directly  bonded  to  it.  This 

cobalt  atom  contributes  0.64  electrons.  However,  the  electron  demand  placed  by  the 

methylene  fragment  is  so  large  that  the  neighboring  atoms  in  the  surface  layer  as 

well  as  metal  atoms  in  the  bulk  also  contribute.  The  metal  atoms  adjacent  to  the 

adsorption  site  contribute  0.19  electrons  each  and  those  in  the  B'  and  C’  layers 

0.32  and  0.20  electrons  respectively.  In  the  methyl  case  charge  transfer  from  the 

inner  metal  layers  is  very  small  (see  21). 

Numerous  examples  of  carbene  complexes  of  type  40  either  substituted  or 

la  11 

unsubstituted  are  known.  ’  Their  bonding  was  analyzed  in  detail  by  our 
group^7a’d  as  well  as  by  others^e  and  we  will  not  dwell  on  this  point  here.  The 
reader  is  also  referred  to  a  recent  relevant  paper  which  reports  the  preparation  art 
the  characterization  of  the  first  simple  unligated  carbene  complex  FeCH2  Note, 
however,  that  to  our  best  knowledge  there  is  not  a  single  known  structure  of  a 
bimetallic  carbene  complex  of  type  40.  These  complexes  exist  commonly  in  a 

bridged  geometry  but  1-fold  structures  have  been  suggested  as  intermediates  in  the:. 
29-31 

reactions . 


We  have  studied  two  bridged  geometries,  perpendicular  56  and  parallel 
: ’ ■  It  is  advantageous  to  analyze  these  geometries  together  although  they 
differ  substantially  in  energy,  56  being  lower  in  energy.  We  will  come  back  to 
this  point  later. 

As  in  the  1-fold  geometry  it  is  useful  to  discuss  the  metal  -  methylene  bond 
in  terms  of  a  and  n  contributions.  The  a  interactions  are  expected  to  be  similar 
the  two  geometries,  as  those  interactions  remain  reasonably  invariant  to  rotation 
around  the  z-axis.  Futhermore,  in  comparison  to  the  on- top  geometry  40,  the 
changes  in  the  a  bonding  are  expected  to  be  similar  to  those  that  occur  in  the 
methyl  systems,  i.e.,  7  -->  8.  For  the  a  framework  we  indeed  find  strong 
similarities  between  the  methyl  and  the  methylene  systems.  In  both  systems  a 
interactions  are  stronger  in  the  bridged  geometry. 


P«rp«ndiculor  PoraiM 
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The  total  density  of  states  for  methylene  on  Co(0001) ,  bringing  in  the 
perpendicular  conformation  56,  is  shown  in  Figure  14.  The  total  DOS  for  the 
parallel  conformation  is  similar,  though  it  differ  in  critical  aspects  to  be 
discussed  below.  Decompositions  of  the  DOS  identify  methylene  levels.  Note  that 
both  a  and  p  bands  are  stablized  in  the  adsorbate  -  surface  complex. 


The  dramatic  difference  between  perpendicular  and  parallel  geometries,  56 
and  57  is  revealed  the  metal-carbon  COOP  curves  which  are  displayed  in  Figures 
15a  and  15b  for  56  and  57  respectively.  Note  that  in  contrast  to  8,  here 
symmetry  dictates  identical  COOP  curves  for  Col-C  and  Co2-C,  of  which  only  one  is 
shown.  The  contribution  of  states  descending  from  the  methylene  p  state  to  M-C 
bonding  are  approximately  twice  as  great  in  56  than  in  57.  In  the 


43 


perpendicular  geometry  the  M-C  bonding  is  roughly  half  n  and  half  a  in  character 
while  in  the  parallel  a  bonding  is  dominating. 

Figure  15  here 

The  stronger  n  bonding  in  56  compared  with  57  results  from  a  more 
efficient  overlap  in  the  perpendicular  geometry  between  the  p  orbital  and  the 
appropriate  metal  states.  In  the  parallel  conformation  the  overlap  between  the 
carbene's  p  orbital  and  the  top  of  the  dz2,  for  example,  is  poor  (see  58)  so 
that  there  is  practically  no  contribution  of  dz2  to  the  p  band.  In  the 
perpendicular  geometry,  on  the  other  hand,  overlap  is  excellent  (see  59)  and  6% 
of  the  dz2  states  contribute  to  the  M-p  bonding  band.  Other  metal  states  that 
contribute  significantly  to  w-bonding  are  the  top  of  the  s-band  (  -5%,  see  60) 
and  the  top  of  the  dx2-^2  band  (  -1%).  dxz  also  contributes,  although  its  share  is 
very  small  (-2%,  see  61).  The  stronger  rr  interraction  in  the  perpendicular 
geometry  are  reflected  in  the  greater  M-C  overlap  populations.  Accompanying  this  is 
a  weakening  of  the  bridged  M-M  bond,  but  the  magnitude  of  that  effect  is  small  (see 
Table  5)  . 


The  stronger  M-C  bonds  in  56  suggest  that  this  conformation  should  be  more 

stable  than  57,  as  observed.  For  all  metal  surface  we  calculate  substantial 

energy  differences  between  the  twc  conformations  (ca.  1 . OeV  for  Co  and  Cr  ca.  0.5eV 

for  Ti).  Note  that  these  energy  differences  also  give  the  barriers  for  rotation  of 

the  CHj  group  around  the  z-axis.  It  is  interesting  to  note  that  these  barriers  are 

similar  in  magnitude  to  those  found  in  molecular  carbene  complexes  of  the  1-fold 
27  29  30 

type.  ’  This  similarity  is  not  surprising  because  the  key  orbital 

interaction  that  dictates  the  magnitude  of  these  barriers  is  similar  in  the 
molecular  and  the  surface  systems. 
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Cur  result  that  the  perpendicular  geometry  is  substantially  more  stable  that. 

the  parallel  geometry  agrees  with  that  of  Minot  ec  al18  but  is  in  contrast  to  the 

conclusions  reached  by  Muetterties  ec  al,  who  suggested  that  the  parallel  geometr 

is  favored  by  strong  me tal -hydrogen  interactions.  As  in  the  previously  discussed 

systems,  we  find  small  antibonding  interactions  between  the  metal  and  the  hydrogen-: 

in  either  56  or  57  and  for  all  electron  counts  (Table  5).  Note  that  our 

results  are  in  line  with  the  known  structures  of  analogous  bridged  molecular 

complexes,  which  invariably  have  a  perpendicular  geometry  with  respect  to  the 

bridged  M-M  vector.  Relevant  molecules  can  be  found  in  references  la  and  11. 

It  is  interesting  to  compare  our  results  with  those  of  Silvestre  and  Hoffmann  who 

have  found  that  the  related  vinylidene  (CH2-C:)  binds  to  a  Pt  surface  most  strong!-.- 

2 1 

in  the  parallel  geometry  62.  The  corresponding  perpendicular  conformation  is 
substantially  higher  in  energy. 
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Although  at  first  glance  this  result  may  seem  to  be  in  conflict  with  that  of 
for  methylene,  in  fact  it  is  not.  Examination  of  both  systems  reveals  that  exact!-.- 
the  same  orbitals  and  interactions  dictate  the  higher  stabilities  of  56  and 
relative  to  the  corresponding  90°-rotated  structures.  The  reversed  final 
geometrical  outcome  results  from  the  different  topologies  of  the  orbitals  of  the  two 
carbenes  at  hand.  In  CH2  the  carbene  acceptor  p  orbital  is  perpendicular  to  the  Hlh 


plane  while  in  vinylidene  it  lies  in  the  molecular  plane.  Thus,  interactions  with 
the  carbene 's  p  orbital  dictate  a  planar  geometry  62  for  the  vinylidene  and  a 
perpendicular  geometry  56  for  methylene.  Note  that  essentially  the  same  reasons 


lead  us  to  construct  the  planar  ethylene  molecule  from  vinylidene  and  hyorogen 
tetrahedral  ( i . e . ,  "perpendicular")  molecule  (methane)  from  methylene  and  H, . 

An  important  result  of  the  calculations  is  that  the  adsorption  of  methylene 
fragment  on  all  metal  surfaces  examined  has  essentially  the  same  energy  in  the 
1-fold  perpendicular  geometry  56  and  in  the  1-fold  on-top  geometry  40.  This 
result  is  to  be  contrasted  with  methyl  adsorption  where  the  1-fold  site  (6)  of  Co. 
Cr  and  Ti  is  preferred  over  the  2-fold  bridged  site  (7)  by  l.leV,  0 . 9eV  and  0.5eV 
respectively.  Why  do  methyl  and  methylene  choose  different  adsorption  sites1  The 
answer  lies  again  in  the  presence  of  the  empty  p  orbital  of  methylene.  We  have 
emphasized  above  for  methyl  that  as  the  organic  fragment  moves  to  a  more  highly 
coordinated  site  its  interaction  with  the  surface  strengthens,  pushing  many  d-stat-c 
to  higher  energies.  The  overall  effect  is  that  the  total  energy  rises  and  the 
bridged  sites  become  less  stable  than  the  1-fold  site.  The  same  effects  are 
operating  in  the  methylene  system  but  here  the  additional  n- stabilization  which  is 
gained  in  the  bridged  geometry  is  so  large  that  it  overrides  the  destabilization 
which  results  from  the  higher  energy  of  the  metal 

centered  bands,  ano  the  total  energy  drops.  Indeed,  a  comparison  of  the  1-fold  sit 

with  the  bridged  parallel  position  where  the  jr-stabilization  is  weak  reveals  the 

’normal"  stability  order,  with  the  1-fold  site  having  the  lower  energy. 

Changing  the  metal  surface  has  a  small  effect  on  the  relative  energies  of  the 

1-fold  and  the  perpendicular  bridged,  sites.  For  Co,  Cr  and  Ti,  the  energy 

differences  are  0.002,  0.11  and  0.02  eV,  respectively,  favoring  56.  Similarly, 

21 

on  a  Pt  surface,  vinylidene  favors  the  2-fold  bridged  site.  These  results  are 

18 

also  in  agreement  with  those  of  Minot  and  co-workers. 

Finally,  we  comment  on  the  charge  distribution.  We  remind  the  reader  that  in 
analyzing  the  charges  it  is  useful  to  look  at  the  separated  molecules  as  [metal, 
and  [organic  fragment]'.  In  the  case  of  methylene  the  charge  -  transferred  extreme 
actually  CHj2',  p  orbital  completely  filled.  As  the  metal  -  adsorbent  interactions 
become  more  efficient  electrons  are  transferred  from  the  adsorbent  to  the  metal. 
The  calculational  results  for  methylene  (see  Table  5)  are  very  similar  to  those 
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obtained  in  Che  methyl  system  (see  34) .  In  both  systems  bridging  increases  the 
interactions  between  the  metal  and  the  organic  fragment  and  consequently  charge  is 
transferred  from  the  adsorbent  back  to  the  metal.  For  example,  the  charge  on  the 
-ethylene  drops  from  -1.54  to  -1.37  when  the  geometry  is  changed  from  the  on- top  r 
the  2-fold  parallel  on  the  Co  surface.  More  interestingly,  the  migration  of  the 
organic  fragment  from  a  1-fold  to  a  2-fold  site  is  accompanied  by  a  massive 
reorganization  of  charge  within  the  metal  slab;  electrons  flow  from  the  surface 
layer  into  the  inner  layers,  avoiding  the  high  lying  bands  which  are  centered  arou 
the  adsorbing  metal  atoms.  It  is  instructive  that  quantitatively  these  inner  meta 
charge  reorganizations  are  practically  identical  for  the  geometry  changes  40  --> 

2Z  (Table  5)  and 

~  -->  8  (see  34).  This  emphasizes  that  similar  basic  changes  occur  in  the 
metal-adsorbent  interactions  when  the  organic  fragment  migrates  from  a  1-fold  to  a 
2 -fold  site,  regardless  of  whether  it  is  methyl  or  methylene.  On  top  of  these  has 
changes  additional  interactions,  specific  to  the  organic  molecule  at  hand,  may  alt 
somewhat  the  final  picture.  Thus  when  the  methylene  is  rotated  to  geometry  56 
and  strong  n - interactions  are  put  into  play,  additional  charge  (relative  to  57) 
flows  from  the  methylene,  which  becomes  less  negative,  into  the  metal  slab.  There 
it  is  partitioned  mainly  between  the  surface  and  the  B'  layers. 


THE  3 -FOLD  BRIDGING  OR  CAPPING  GEOMETRY.  42 

The  analysis  of  the  bonding  at  this  site  is  straightforward,  using  the 
experience  gained  in  the  study  of  the  methyl  system  and  of  the  various  adsorption 
sites  of  methylene.  The  COOP  curves  of  the  Co-C  bonds  are  displayed  in  Figure  16 

Figure  16  here 

The  numbering  of  the  metal  atoms  and  the  precise  orientation  of  the  methylene  grc 
is  shown  in  63.  The  solid  line  in  Figure  16  gives  the  COOP  curve  for  the  Col-C 
bond  and  the  dotted  line  shows  the  corresponding  curve  for  the  two  identical  Co2- 
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Figure  16:  COOP  curves  for  the  Ml-C  and  M2-C  Bonds.  CH2  group  Is  above  the 
triangular  hollow  of  the  Co(0001)  surface. 


and  Co3 -C  bonds.  The  COOP  curve  shows  chat,  similarly  to  the  2 -fold  56  system, 
also  in  63  both  a  (band  at  -13.6eV)  and  n  interactions  (band  at  -12.2eV) 
contribute  effectively  to  Co-C  bonding.  Typical  a  and  it- interactions  with  the 
ratal's  d,2  band  are  depicted  in  64  and  65  respectively.  Similar 

interactions  with  other  metal  bands  in  particular  s  and  to  some  extent  dxz  and  d.._ 
also  contribute  to  bonding. 
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Figure  16  reveals  an  interesting  aspect  of  the  bonding  at  the  3-fold  site  and 
provides  a  nice  demonstration  of  the  analytical  power  of  COOP  curves  in  general, 
j-bonding  (peak  at  -13.6eV)  is  contributed  equally  (except  for  a  small  perturbatio: 
induced  by  the  non- symmetrical  arrangement  of  the  methylene  hydrogens)  by  all  thre. 
metal  atoms.  In  contrast  Co2  and  Co3  contribute  to  ^-bonding  (peak  at  -12.2eV)  mu 
more  than  Col.  This  result  can  be  interpreted  quite  easily.  In  64  the  p 
orbital  of  methylene  Is  aligned  parallel  to  the  Co2-Co3  vector  allowing  strong 
overlap  with  the  bands  centered  around  these  metal  atoms.  The  topology  of  the 
interacting  orbitals  is  similar  to  that  in  the  perpendicular  2-fold  geometry.  The 
interaction  of  Col  with  the  p  orbital  is  much  smaller,  the  situation  resembles  tha 
of  Col  in  the  parallel  2-fold  geometry.  Upon  rotation  of  the  methylene  group  by  3 
in  a  clockwise  direction  Col  and  Co3  change  roles;  now  Col  and  Co2  contribute  to 
T-bonding  more  than  Co3 .  At  intermediate  rotation  angles  metal  bands  combine 
effectively  to  maintain  the  level  of  w-bonding,  so  that  at  the  triangular  hollow 
site  the  adsorbed  methylene  fragment  is  essentially  freely  rotating. 

We  have  already  indicated  in  various  places  along  the  discussion  that  when  th 
organic  fragment  migrates  to  a  more  highly  bridged  site  its  interactions  with  the 
metal  strengthen.  Geometry  42  is  no  exception.  For  all  three  metal  surfaces 
the  total  M-C  overlap  population  (Table  5)  is  higher  in  the  3-fold  geometry  (eg 


tor  than  m  the  2-fold  perpendicular  geometry  (0.792  for  Co )  .  Other 

trends  that  we  have  observed  for  the  8  - ->  9  migration  of  methyl  hold  also  for 
T.ethvler.e  (Table  A).  As  most  of  the  energy  which  is  associated  with  n  bonding  was 
alreadv  gained  in  the  2-fold  geometry  further  bridging  leads  to  the  "normal" 
increase  in  the  total  energy,  or  equivalently  to  a  decrease  in  the  binding  energy  t 
the  organic  molecule  to  the  metal.  Binding  of  methylene  at  the  3- fold  site  of 
cobalt  is  consequently  by  0.6eV  less  favorable  than  binding  at  the  2 -fold  or  at  the 
on-top  sites.  Similar  trend  holds  for  the  Cr  and  the  Ti  surfaces. 

Finallv,  we  note  that  we  are  not  aware  of  molecular  trinuclear  transition  met. 


carbene  complexes  that  are  stable  in  the  3- fold  capping  geometry.  Polymetallic 


carbene  complexes  usually  adopt  the  2-fold  perpendicular  geometry. 


la,  11 


We  can  now  proceed  to  analyze  the  adsorption  of  the  last  fragment  -  methvr.e 


ADSORBED  METHYNE.  CH 

We  will  compare  three  different  adsorption  sites  66,  67  and  68. 

Again,  the  calculations  are  performed  for  an  "unbiased"  C-M  distance  of  2 . 1&  for 
all  three  geometries.  The  binding  energies,  overlap  populations  and  charges  are 
tabulated  in  Table  7.  We  will  analyze  the  bonding  based  on  DOS  and  COOP  curves. 


The  MO’s  of  the  methyne  group  before  and  after  the  adsorption  on  the  hexagor.a 
metal  surface  in  a  capping  mode  are  shown  in  Figure  17.  The  reason  we  look  at  the 
capping  geometry  is  that  this  site  is  the  preferred  one  for  high  metal  electron 
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counts.  Ac  the  left  side  are  Che  surface  scares  plus  che  CH  levels  before  the 
chemisorpcion  cakes  place,  Che  righc  side  is  whac  happens  afcer  CH  is  chemisorbed 

Figure  17  here 

che  me  Cal  surface.  CH  has  five  MO's:  one  of  Che  Cwo  carbon  sp  hybrids  pointing 
coward  che  hydrogen  forms  C-H  a  and  a*  orbicals.  The  ocher  hybrid,  direcced  awav 
from  che  hydrogen,  is  mainly  non-bonding,  and  so  are  Che  degenerace  p  p.,  orbital 
The  a  orbical  is  low  in  energy  and  we  expecc  it  interaccs  poorly  wich  the  rr.ecal  d 
band.  a*  is  high  up  in  energy  (off-scale  in  Figure  17),  also  out  of  effective 
interaction  range.  What  is  left  for  bonding  with  the  surface  are  the  n  ar.d  p  ,  r... 
orbitals . 

Lee  us  focus  on  Figure  17b.  We  see  that  n  is  pushed  down  by  -leV;  at  higher 
energy  p  is  smeared  out  over  the  region  of  -12.3  -  -8eV.  Both  CH  orbitals  interac 
with  che  metal  surface.  Which  metal  orbitals  are  effective  in  this  interaction? 
Figure  18  shows  projections  of  metal  dz2  and  dxz ,  d^z .  The  dz2  orbical  picks  up  a 


resonance  wich  the  CH  p  orbital  at  -12eV:  interactions  of  type  69a  are 
involved.  The  small  peak  in  the  p  projection  in  Figure  17  at  -8eV  is  also  a 
consequence  of  this  mixing.  The  dxz ,  d^z  set  has  resonances  with  both  CH  p  and  r. 
Figure  18.  Combinations  of  type  69b  and  69c  are  effective. 


What  happens  as  CH  moves  to  other  adsorption  sites,  on- top  or  bridging?  We 
r.ot  show  che  DOS  plots  here,  but  the  n  level  is  pushed  down  a  lictle  less  chan  in 
che  capping  sice.  More  of  a  difference  is  observed  in  the  fate  of  the  p  level.  : 
is  at  -  -12.5eV  for  68,  -12.1eV  for  67  and  -11.7eV  for  the  on- top  geometrv. 

Motivated  by  Muetterties  and  co-worker's  speculation  that  CH  group  might  r.ot 
perpendicular  to  the  surface  due  to  M-H  interaction  ,  we  compared  two  geometries 
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Figure  17:  DOS  of  the  CH  +  Co(0001)  slab  system  before  (a)  and  after  (b)  CH 
chemisorbs  on  the  metal  surface  (capping  geometry) .  CH  states  are  shown  by 
horizontal  bars  (a)  or  darkened  area  (b) . 


’2^  "  67  anc*  ZZb-  Our  calculation  shows  no  significant  energetic 
difference,  2Qb  being  slightly  higher  in  energy.  Nearly  all  the 

organometallic  analogues  have  CH  or  CR  in  a  symmetrical,  or  capping  geometrv. 
21. Binuclear  complex  with  a  bridging  CR  are  well-know*.,  as  are 
mononuclear  L  M"CR  acetylene  analogues . 3 ' ^ 


a  70  b  71 

To  summarize  our  results  concerning  various  CHx  groups  on  the  surface,  we  cai 
simply  state  that  for  equal  C-M  distances,  for  all  adsorption  sites,  the  CHX  grou; 
would  prefer  to  adopt  a  geometry  which  can  restore  it's  missing  C-H  bond(s).  Thu 
CH^  on-top,  CH9  bridging,  and  CH  capping  are  energetically  the  most  stable 
geometries.  This  agrees  with  previous  theoretical  studies  and  the  available 
experimental  information. 

So  far  we  have  said  nothing  of  reactivity.  But  most  of  the  experimental 
studies  of  CHX  groups  cannot  avoid  the  dynamic  aspect.  Reactions  occur.  For 
example,  CH2  on  Fe(100)  is  believed  to  dissociate  into  CH  +  H2 ,  beginning  at  the 
bridging  position,  with  one  of  the  C-H  bonds  parallel  to  the  surface  j3  On  a  Co 
surface  the  dehydrogenation  of  CH2  occurs  at  a  relatively  low  temperature.  180K, 
indicating  a  small  barrier  for  the  process  CH2acjs  - ->  CHacJs  +  Hads  ^  t^‘e  nex 


section  we  are  going  to  analyze  the  migration  of  CHX  groups  on  a  surface  and  how 
electronic  factor  affects  various  reaction  processes. 


Migration  on  the  Surface 

Migration  of  species  adsorbed  on  a  surface  is  a  phenomenon  of  substantial 
significance  in  many  areas  of  surface  science  (e.g.,  adhesion,  lubrication,  etc.;. 
For  catalysis,  the  case  of  interest  here,  knowledge  of  the  barriers  for  the 
migration  of  the  adsorbed  fragments  is  fundamental  to  the  understanding  of  the 
mechanism  of  surface  reactions.  In  the  FT  synthesis,  for  example,  it  is  important 
to  know  if  in  the  coupling  step  (step  3,  scheme  1)  the  rate  determining  process  is 
the  migration  of  the  two  (or  more)  fragments  towards  each  other,  or  if  it  is  the 
chemical  barrier  associated  with  the  coupling  reaction  itself.  Unfortunately,  the 
experimental  characterization  of  surface  migration  (or  diffusion)  is  difficult,33 
and  for  the  CHX  fragments  little  such  information  is  available. 

Recently,  Muetterties,  Shustorovich  and  Baetzold  developed  a  remarkably  useful 
and  versatile  theoretical  model  for  surface  migration  of  atoms  and  diatomic 
molecules.  This  is  based  on  a  Morse  potential  for  the  metal  -  adsorbent  interactions 
and  on  the  assumption  that  the  M-X  (M  -  metal,  X  -  adsorbent)  bond  order  is 
conserved  along  the  migration  path.36  Our  approach  is  different  and  is  based  on  th 
actual  calculation  of  several  selected  points  on  the  potential  energy  surface.  All 
the  necessary  data  for  analyzing  the  migration  of  the  CHX  fragments  on  the  metal 
surface  have  been  already  presented  above.  A  disadvantage  of  our  approach  has 
already  been  mentioned  above  -  we  rely  on  a  molecular  orbital  method  which  is  very 
approximate,  especially  when  it  deals  with  the  energetics  of  bond  formation.  The 
conclusions  we  reach  gain  a  little  strength  when  they  are  supported  by  overlap 
population  analysis,  which  together  with  the  energies  and  ocher  bonding  arguments 
form  a  self-consistent  conceptual  picture. 

In  drawings  £2 ,  73  and  74  we  have  collected  what  the  computations 
give  for  the  relative  binding  energies  (in  eV)  at  the  three  symmetry  distinguished 
adsorption  sites,  of  methyl,  methylene  and  methyne ,  respectively  In  each  drawing 
we  present  the  calculated  values  for  Co,  Cr  and  Ti  surfaces 
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A  low  migration  barrier  can  be  observed  only  if  there  are  at  least  two  sites 
connected  or  nearby  on  the  metal  surface  for  which  the  binding  energies  of  the 
organic  fragment  are  similar.  Let  us  call  these  sites  A  and  B.  In  this  case  the 
lowest  energy  path  for  the  migration  of  the  fragment  is  simply  along  the 
A  -->  B  -->  A  ...  path.  The  energy  barrier  for  the  migration  is  then  determined  by 
the  energy  differences  between  points  A  and  B.  The  relative  energies  of  the  other 
binding  sites  (C,  D,  etc  )  are  of  no  significance  to  the  migration  process  because 


even  if  these  points  are  very  high  in  energy  they  can  be  circumvented  bv  migration 
along  the  A  -->  S  -->  A  path. 

Let  us  look  first  at  the  migration  of  a  methyl  fragment  that  is  shown  in 
22.  For  a  cobalt  surface  the  site  with  the  highest  binding  energy  (i.e  the 
lowest  point  on  the  potential  energy  surface)  is  the  1-fold  "on-top"  site  The  two 
other  sites  are  considerably  higher  in  energy.  The  lower  energy  one  of  these  is  th 
2-fold,  l.leV  higher  in  energy  than  the  1-fold  site.  The  theoretical  prediction  is 
that  the  migration  of  methyl  on  a  cobalt  surface  will  require  relatively  high 
temperatures  and  should  proceed  mainly  via  the  path  shown  by  the  dotted  line  in 
21-  Migration  occurs  along  the  (formal)  metal-metal  bonds.  When  the  electron 
count  on  the  metal  is  lowered  the  bridged  positions  are  differentially  stabilized 
and  the  barriers  to  migration  should  be  considerably  lower.  For  the  Ti  surface  the 
barrier  to  migration  along  the  path  shown  in  76  (which  is  the  lowest  in  energy) 
is  only  O.leV.  For  the  Cr  surface  the  barrier  is  calculated  to  be  the  same  <0  9eV) 
along  either  path  75  or  76. 


liras 

am 


Why  does  the  migration  barrier  decrease  as  the  metal  d  band  is  depopulated? 
Remember  the  main  interaction  is  between  the  methyl  lone  pair  n  and  the  metal  d 
band.  77.  The  metal  d  band  region  is  weakly  C-M  antibonding  (c.f.  Figure  6). 

The  strength  of  the  antibonding  character  depends  on  the  overlap  between  n  and  the 
metal  d  band  The  higher  the  coordination  number  of  the  CH^  group,  the  stronger  the 
C-M  antibonding  character  in  the  metal  d  band.  Thus  at  high  d  band  filling  the 
capping  geometry  is  less  stable,  because  of  its  stronger  C-M  antibonding 
contribution  within  the  increasingly  occupied  d  band.  At  lower  band  filling  the  C-M 


antibonding  feature  is  reduced,  and  the  higher  coordination  sites  of  CH-,  stabilize  n 


more.  So  the  barrier  is  reduced 


T 


For  methylene  the  situation  is  different.  Migration  on  a  cobalt  surface  occurs 
with  no  barrier  along  a  path  (with  no  rotation  of  the  CHj  group)  shown  by  the  dashec 
line  in  78.  The  same  lowest  energy  migration  path  is  calculated  for  CH0  on  the 
Cr  and  Ti  surfaces  (see  also  the  more  detailed  discussion  above). 


The  migration  of  a  CH  fragment  is  relatively  facile;  the  barriers  along  path 
'9  are  only  0.3eV,  0.4eV  and  O.leV  for  Co,  Cr  and  Ti  respectively.  This 
migration  path  remains  the  lowest  in  energy  for  all  metals,  and  for  those  with  lower 
electron  counts  the  binding  preference  for  the  2-  and  3-fold  sites  even  increases 


The  reason  that  the  migration  barriers  for  CH2  and  CH  do  not  change  much  at  lov 
metal  band  filling  derives  from  the  fact  that  in  addition  to  orbital  n  (a)  some 
other  orbitals  (p)  also  play  a  decisive  role  in  the  bonding.  Details  of  the 
analysis  have  been  given  above.  One  can  say  that  the  p  orbital,  being  very  close  t: 
the  d  band  in  energy,  mixes  very  strongly  with  the  metal  orbitals  (c.f.  Figures  12b. 


15 , 19  and  21)  . 


This  mixing,  shown  in  80,  may  contribute  to  C-M  bonding  even  at 
high  d  band  filling.  Thus  two  types  of  interaction  77,  80  compete  to 
determine  the  migration  barrier. 


80 

Knowledge  of  the  barriers  to  migration  has,  of  course,  interesting  implications 
regarding  the  coupling  steps  in  the  FT  synthesis.  For  example,  the  calculations 
suggest  that  on  a  cobalt  surface,  the  coupling  of  methyl  and  methylene  may  occur 
preferably  by  migration  of  CH2  towards  the  adsorbed  methyl  rather  than  by  migration 
of  methyl.  Furthermore,  the  coupling  of  two  methyl  fragments  may  involve  a 
substantial  barrier  for  migration,  while  the  coupling  of  two  methylenes  should 
require  no  such  barrier.  To  demonstrate  the  implications  of  these  results  let  us 
assume  that  the  barriers  to  the  coupling  processes  themselves  ar~  small  so  that  the 
migration  processes  determine  the  rate  of  the  FT  reaction.  The  above  results 
suggest  that  in  such  a  case  the  rate  of  coupling  of  two  methyls,  but  not  that  of  two 
methylenes,  should  be  reduced  at  lower  temperatures.  Thus,  the  coupling  of  two 
methylenes  (or  of  a  methylene  and  a  methyl)  would  be  favored  at  lower  temperatures 
and  higher  yields  of  ethylene  expected.  Similar  analysis  can  lead  to  interesting 
conclusions  regarding  the  Cr  and  Ti  metal  surfaces  and  the  other  coupling  reactions. 


The  coupling  reactions  in  Scheme  3  determine  Length  and  olefin  proportion  of 


the  FT  product.  Do  electronic  factors  control  these  reaction?  Why  does  the 

Schuit - Florv  distribution  change  when  a  different  metal  catalyst  is  used?ic  Such 

information  could  certainly  help  chemists  design  new  catalysts  with  specified 
52 

selectivities .  We  do  not  intend  and  certainly  cannot  answer  all  those  questions, 
but  we  hope  to  extract  some  useful  information  from  our  calculations .  In  this 
treatment  we  assume  that  methyl  and  methylene  are  reasonaly  good  models  for  other 
adsorbed  alkyl  (CnH2n.1)  or  alkylidene  (CnH2n)  fragments. 
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Scheme  3 


CH-,  +  CH-, 

All  the  previous  calculations  were  performed  assuming  a  1/3  coverage  for 
reasons  of  computational  economy.  If  we  bring  another  CH^  group  into  the  unit  cell 
the  high  2/3  coverage  will  certainly  lead  to  great  crowding.  A  geometry  such  as 
31  is  forced.  In  our  calculation  for  this  and  all  other  reactions  we  delete  the 
interaction  of  the  groups  between  unit  cells,  so  that  there  is  no  crowding  whose 
cause  is  inter-cell  interaction.  In  other  words,  we  are  modeling  a  reaction  at  ver 
low  coverage. 
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The  DOS  curves  for  CHj  on  Co  in  such  a  coverage  (not  shown  here)  show  a 
splitting  of  all  CH2  states.  Such  a  splitting  in  the  occupied  states,  shown  ir 
32  for  a  lovlving  a  orbital,  is  highly  destabilizing.  This  is  a  typical  four- 
electron  two-center  repulsion.  The  COOP  curves  also  indicate  regions  of  C-C 


0^-0  ch^>-o 


antibonding.  83  gives  computed  binding  energies  per  CH2 ,  down  considerably  from 


the  low  coverage  case  (Table  4) .  The  differences  in  the  binding  energies  in 

83  and  Table  4,  which  we  call  the  proximity  barrier,  are  0.77eV,  0.53eV,  and 


-0.05eV  for  Co,  Cr  and  Ti. 
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What  happens  when  two  CHj  groups  couple?  The  reaction  begins  withe  both  CH2 
lone  pairs  nearly  filled,  i.e.  near  a  CH2 ’  representation.  A  new  C-C  a  bond  forms. 

if 

and  as  usual  we  must  consider  a  and  a  combinations,  n^±n£ ■  Both  are  filled 
initially,  but  as  the  C-C  bond  is  more  completely  formed  the  a *  combination  will  be 
pushed  up  and  eventually  will  dump  its  electrons  into  the  metal  d-band. 

Let  us  sample  such  a  reaction  on  the  surface  by  using  a  transit  along  a  single 


reaction  coordinate  8  84,  The  "umbrella  handle"  of  the  CH2  group  is  moved  along 
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ur.d  perpendicular  to  an  arc  of  a  circle.  The  circle  is  defined  as  "standing''  or.  - 
petal  atom,  such  that  when  the  rotation  angle  9  -  90°  the  C-C  distance  is  1.5-  3 
i C  to  surface  separation  is  2.59&).  9  -  0°  corresponds  to  the  initial  geometry 

with  the  "umbrella  handle"  perpendicular  to  the  surface  and  bonded  to  the  metal 
atom.  This  transit  may  not  be  optimal,  a  true  reaction  path,  but  it  can  give  us 
some  information  concerning  the  controlling  electronic  factors.  Full  exploration  of 
the  potential  surface  is  not  realistic  with  our  computational  capabilities. 

Figures  19  and  20  show  the  methyl  lone  pair  contribution  to  the  total  DOS  and 
C-C  COOP  curve  along  the  reaction  coordinate  9.  Before  the  coupling  starts  (9  -  0°) 
n-CHj  already  has  split  into  two  peaks  around  -13eV.  The  higher  one  is  n^-n2  (c.f. 
82)  and  the  lower  one  n^+n2 .  But  at  a  C-C  separation  of  2.51&  n]/n2  i-s  still 
hardly  antibonding.  At  9  -  30°  the  antibonding  portion  (-12  -->  7eV,  of  Figure  19) 
separates  from  the  bonding  part  (-  -14eV)  and  at  9  -  60°  is  well  above  the  Fermi 
level . 

Figures  19,  20  here 

What  if  the  position  of  the  Fermi  level  changes?  From  the  above  discussion  the 
total  energy  of  the  system  increases  along  the  path  as  n^-n,,  becomes  more 
antibonding.  When  n^-n2  is  pushed  above  the  Fermi  level  it  becomes  empty.  After 
this  turning  point  in  the  reaction  path  the  total  energy  should  go  down  again  as 
n^+n2  becomes  more  bonding.  The  position  of  the  Fermi  level  determines  the  turning 
point  and  the  barrier  height.  The  higher  the  Fermi  level,  the  further  away  that 
turning  point  from  the  starting  geometry,  and  the  higher  the  barrier.  Figure  21, 
showing  the  computed  energy  profile  along  the  idealized  reaction  path,  confirms  this 
line  of  reasoning.  The  coupling  barrier  on  Co  is  lower  than  that  on  Cr,  which  in 
turn  is  lower  than  that  on  Ti. 

Figure  21  here 


£H3  ±  (?H2 

This  is  the  chain  propagating  step  in  the  FT  process.  The  binding  energies  E1 


in  85  tell  us  that  some  55  kcal/mol  per  CH-j  +  CH2  (the  proximity  barrier)  is 


ction  c 


needed  to  Siring  the  C.  group  :o  the  geometry  S3,  compared  to  36  Kcal  mol  0. 
per  CH^)  for  S3. 
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The  reaction  coordinate  is  chosen  again  to  be  a  single  parameter  9,  86. 

The  C-C  distance  (decreasing  from  2.51  to  1.548),  the  angle  a  between  the  C-C  bond 
and  the  CH^  "umbrella  handle"  (also  decreasing,  from  90°  to  0°)  and  the  HCM2  angle  3 
(from  120°  to  109.5°)  on  the  CH2  group  are  chosen  to  vary  linearly  with  9.  9  itself 

changes  from  90°  to  109.5°. 


2.31  A 

1 

j  u 

Xs  If1' 

86 

Along  this  transit  we  expect  that  one  C-M  bond  (from  CH^  to  a  metal  atom)  is 
broken  and  C-C  bond  formed.  The  CH^  group,  formally  bonded  to  the  surface  by  a 
double  bond,  loses  some  of  its  C-M  bonding  when  the  hydrogens  are  bent  downwards. 
This  is  because  the  CH2  orbitals  rehybridize  and  in  the  final  geometry  there  is  a 
hybrid  pointing  away  from  the  metal  atom,  shown  in  87.  That  hybrid  is  gradually 
pushed  up  by  an  approaching  methyl  n  orbital,  88.  Thus  the  potential  energy 
curve  should  rise  due  to  the  four-electron  repulsion  88,  and  then  drop  down 
after  the  antibonding  level  is  pushed  above  the  Fermi  level  and  transfers  its 


electrons  to  the  metal. 


Figure  22  shows  the  contribution  of  the  methylene  p  orbital  (in  local 
coordinates')  to  the  total  DOS  along  the  reaction  coordinate.  A  portion  of  the  tota 
contribution,  corresponding  to  the  C-C  a  level/ gradually  climbs  from  below  the 
Fermi  level.  Figure  23  is  the  computed  potential  energy  curve  along  the  reaction 
path.  Again  as  the  metal  Fermi  level  is  lowered  the  barrier  height  is  reduced  and 
the  peak  shifted  to  the  left  side,  for  reasons  similar  to  those  discussed  for  the 
CHj+CHj  case.  The  n-CH^  and  n-CH2  states  also  have  a  portion  corresponding  to  the 

■fa 

C-C  a  level  rising  up  along  the  reaction  path  as  expected,  but  these  graphs  are  no 
shown  here . 

Figures  22 .  23  here 


ch2_l_ch2 

The  surface  precursor  could  start  in  either  a  parallel,  89,  or  a 
perpendicular  geometry,  90.  The  possibility  of  the  perpendicular  mode  is 
anticipated  from  our  earlier  discussion  that  the  capping  geometry  for  a  CH2  group  o 
the  hexagonal  surface  may  have  lower  energy  for  low  d  band  filling  (see  78). 

Let  us  see  how  the  calculations  check  out  our  expectations. 


Figure  22:  Methylene  p  state  evolution  along  a  reaction  path  on  Co(0001). 
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Figure  24  shows  Che  CH0  contribution  to  the  total  DOS  for  geometry  89a. 

The  calculation  is  performed  assuming  a  low  coverage,  i.e.,  interaction  between  Ci 
groups  other  than  those  two  under  consideration  is  deleted.  The  two  CH0  groups  a' 
C-C  separation  of  2.51&  interact  with  each  other  weakly,  so  each  MO  level  is  spli' 
into  two.  as  we  can  see  from  Figure  24.  The  "four-electron  repulsion"  reduces 
somewhat  the  stability  of  the  system.  91  shows  the  binding  energies  for  each 
CH2 .  The  numerical  values  tell  us  that  it  costs  some  12  kcal/mol  energy  (proximi 
barrier)  to  bring  two  CH2  groups  together  at  a  separation  of  2.5lX. 


V  V 
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E  per 

CHj  (»v: 

e, 

Co 
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004 

Cr 
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Ti 
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91 

Figure  24 

hero 

The  other  geometry,  90a,  wi-~h  each  CH2  group  sitting  directly  above  the 
metal  triangular  hollow,  requires  a  C-C  separation  of  1.458.  This  is  obviously 
not  a  realistic  starting  geometry  and  our  calculation  indeed  gives  a  large 


Figure  24:  Total  DOS  (d 
CH2  +  CH2  on  Co(0001)  sy 


■.e  ar.s  . 


i  'r.  To  :  ir.d  the  appropriate  reaction  path  for  that  geometry  is  be  von  i  ,  ;r 
instead  we  will  explore  some  limited  parts  of  the  potential  surface,  in 
articular  those  verv  close  to  the  final  product  in  the  coupling  reaction. 

cut  what  is  the  final  geometry?  As  we  have  said  before  the  molecular  corr.cle: 


geometry.  The  C-C  distance  is  -  1.528,  character istic  of 


i  smg.e  tone. 


The  C-M  distance  is  -  2.0.8  (M-Fe)^a,  s 


uggestmg  a  single 


’  t*>  "  1  a  I*'  .  M 


von.d.  And  the  hydrogens  bend  awav  from  the  metal. 


Vv_Lc-> 

/  \ 

W - M 


Some  information  is  available  concerning  the  geometry  of  chemisorbed  ethvlen 

Ibach  and  Lehwald  have  compared  the  vibrational  frequencies  of  02^  on  Pt(lll)  wi 

those  of  Zeise's  salt  (  K '  PtC  1  ^  '  H2O )  and  concluded  that  the  carbon  atoms  hav 

sp  -hybridization.  Comparison  of  experimental  photoemission  data  with  an  SCF-L 

calculation  on  ethylene  in  different  geometries  by  Demuth  suggested  that  the  C-C 

distance  is  1 . 34  -  1.498  and  the  C-C-H  angle  is  106°  -  109.5°  for  ethylene  on 

Cuf  111)  .  N’i(lll),  Pd(lll),  and  Pt(lll).^  A  recent  HREELS  study  by  Stroscio,  Hare 

3  4 1 3, 

do  also  suggests  that  adsorbed  ethylene  is  characterized  by  sp  hybridization, 
iom.parison  of  NEXAFS  of  gaseous  and  chemisorbed  (on  Pt(lll))  ethylene  led  to  a 
proposal  of  a  single  C-C  distance  of  1.538.^^ 

interestingly,  on  Pt(lll)  precovered  with  oxygen  there  is  a  mixture  of  both 
di -a  bonded  (sp^)  and  m  bonded  (sp^)  ethylene.4^3  On  Ni0(100)  evidence  exists  fo 
two  adsorption  states  of  ethylene.  At  high  temperature  ethylene  also  transform 

into  ethylidyne  H^C -C- . ^  ^  Other  possibilities  include  hvdrogenolvs is44  or 
dehydrogenation  v  of  ethylene  on  the  surface.  Since  the  nioneering  work  of  Rose 
and  Rhodin,14^^  Anderson, 4^e  and  Demuth,4^  there  have  been  many  theoretical  stud;, 
of  ethylene  interacting  with  metallic  clusters  of  finite  size . 4  3  A  re.er.t  ; 


.  .  4  5  f 


studies  ethylene  cn  (111)  layers. 


From  the  available  experimental  information  we  choose  the  : 


93:  C-C  =  1.45A,  the  HCH  plane  bends  upward  by  30°  but  the  HCH  angle  remains 


120°.  We  are  going  to  compare  three  different  adsorption  sites  94,  95,  and 

96,  assuming  a  C-M  distance  2.1&  for  all.  Since  most  experimental  studies  are 

43r 

on  late  transition  metals  ,  we  will  perform  the  calculation  only  for  Co (0001). 
The  choice  of  these  geometries  is  based  on  our  knowledge  of  organometallic 

47 

compounds,  in  particular  those  in  ref.  38  and  the  structure  of  Zeise's  salt. 
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Figure  25  here 


For  the  small  degree  of  bending  or  puckering  assumed,  the  ethylene  orbicals 
will  change  little  on  going  from  the  planar  molecule  to  the  puckered  93.  Figure 
25  confirms  this  and  also  serves  to  remind  us  of  the  orbitals  of  ethylene,  easily 
related  to  those  of  two  interacting  CH2  groups.  ir  and  w  ,  responsible  for  the 
important  forward  and  back-donation  in  the  Dewar-Chatt-Duncanson  Model, ^  acquire 
little  a  character.  Another  a  orbital,  with  its  lobes  directed  inward,  may  also 
play  an  important  role  in  the  bonding,  especially  in  geometry  95.  97  might 

be  an  effective  overlap  in  that  geometry. 


We  compare  the  ethylene  orbital  contribution  to  the  total  DOS  for  those  three 


absorption  sites,  94,  95  and  96  in  Figure  26.  As  we  can  see  upon 
adsorption  all  occupied  orbitals  (up  to  jt)  remain  approximately  where  they  were 
before  adsorption.  The  unoccupied  orbitals  (above  x  )  are  all  smeared  out,  more  so 
for  geometry  96. 

Figure _ 26  here 

Let  us  zoom  in  at  those  three  orbitals  a,  x  and  x  whose  lobes  point  toward  the 
metal  surface  and  thus  should  have  the  greatest  potential  for  metal -ethylene 
interaction.  The  calculations  (summarized  in  part  in  Table  8)  show  that  a  interacts 
little,  most  in  the  capping  geometry  96.  Even  there  it  donates  no  more  than  4% 
of  its  electron  density  to  the  surface.  A  typical  interaction  involved  is 

Table  8  here  | 

98.  A  similar  conclusion  holds  for  ir',  in  the  capping  geometry  x  donates  more 
electrons  and  is  pushed  down  more.  But  in  the  butterfly  geometry  95  there  is  no 
indication  that  the  interaction  is  stronger  than  in  the  sawhorse  geometry,  mostly 
due  to  the  fact  that  x  has  its  lobes  pointing  outward,  away  from  the  metal  atom, 

99.  Table  8  shows  that  x  has  donated  substantial  density  to  the  surface. 


98 


99 


The  biggest  difference  occurs  in  the  x*  orbital,  and  for  this  one  we  offer  a 
decomposition  of  the  DOS  in  Figure  27.  x*  is  smeared  out  more  along  the  energy 
scale  as  the  geometry  goes  from  94  (sawhorse)  to  95  (butterfly)  to  96 
(capping) .  In  the  sawhorse  geometry  x*  interacts  only  with  the  top  of  the  d  band 
(M-M  antibonding).  A  characteristic  interaction  ( dz 2 )  is  shown  in  100.  In 
the  butterfly  modes  it  interacts  with  the  entire  dxz  band  101,  especially 
those  portions  that  match  the  energy  of  x  .  So  we  expect  the  interaction  to  be 


scronger  and  x*  spread  out  more.  The  capping  geometry  96,  still  remains  the 
most  effective  for  interaction  and  we  see  x  dispersed  most.  Again  this  is  due  to 
the  fact  that  more  metal  atoms  and  orbital  combinations  are  available  to  interact 
with  x* .  From  its  band  width  and  proximity  to  the  Fermi  level  we  can  conclude 
that  *T*  has  the  strongest  interaction  with  the  metal  and  is  primarily  respon¬ 
sible  for  the  binding. 


The  top  of  the  wide  x  band  is  C-M  antibonding  and  the  bottom  C-M  bonding.  The 
stronger  the  interaction  between  x*  and  metal  the  more  profound  the  antibonding  or 
bonding  character.  Thus  at  low  band  filling  that  geometry  which  allows  stronger 
interaction  has  more  bonding,  but  at  high  band  filling  more  antibonding.  This  kind 
of  reasoning  can  extend  our  considerations  to  surfaces  other  than  the  one  considered 
here . 


Table  8  collects  the  bonding  information  for  all  three  geometries.  As  we  have 
said  before,  the  x*  interaction  with  the  metal  strengthens  on  going  from  sawhorse 
(94)  to  butterfly  (95)  to  capping  (96)  geometries.  The  stronger 
interaction  for  the  capping  geometry  is  supported  also  by  the  fact  that  ethylene 
donates  more  changes  to  the  metal  (see  last  column  in  Table  8).  Our  calculation 

it 

also  suggests  that  on  the  sawhorse  geometry  x  is  filled  more  and  ethylene 
dissociation  into  Cl^’s  should  be  more  facile. 

The  sawhorse  mode  is  categorized  as  the  di-o-bonding  and  the  butterfly  as  the 
x-bonding  mode  in  surface  science  literature.  Experimental  studies4"*15  as  well  as 
theoretical  calculations4^  tell  us  that  x  is  pushed  down  more  in  the  x -bonding  than 
in  the  di-<7-bonding  geometry.  Our  calculation  is  no  exception,  but  in  addition 
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indicates  the  driving  force  behind  the  stronger  interaction  in  the  butterfly 
geometry . 


CONCLUSIONS 

Let  us  summarize  several  important  results  and  discuss  the  possible 
consequences.  In  the  C^  fragment  part  we  concluded  that  n-CH-j  or  a-C H2 
(-  n-Ct^  )  is  responsible  to  the  binding  of  the  corresponding  species  to  the 
surface.  n-CH^  or  a-Cl^  is  pushed  down  in  energy  and  its  bonding  partner,  the  metal 
d  band  region,  is  characterized  by  weak  M-C  antibonding.  This  is  shown 
schematically  in  103 . 


102 


We  saw  what  happens  as  the  methyl  group  is  moved  from  the  on- top  site  103 

i 

to  the  2 -fold  bridging  site  104  to  the  3 -fold  bridging  site  105.  For 

symmetry  and  overlap  reasons  the  antibonding  feature  in  the  d  band  region  becomes 

stronger  and  the  total  energy  of  the  system  goes  up  along  with  this  migration. 
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Also  Che  methyl  Co  metal  antibonding  character  grows  with  the  metal  d  band  filling. 
The  binding  energy  decreases  for  most  of  the  binding  sites  and  at  the  same  time  the 
difference  between  the  binding  energies  for  these  absorption  sites  is  enhanced  (cf. 
Table  4  and  79) .  Therefore  the  methyl  group  should  become  less  mobile  on  the 
surface  as  d  bands  are  more  occupied.  For  the  same  reason  the  "covalent" 
contribution  E2  to  the  binding  energy  of  CH2  and  CH  on  a  metal  surface  also 
decreases  as  the  metal  d  band  is  more  populated.  The  reduced  mobility  as  the  metal 
d  band  is  more  occupied  contributes  to  a  smaller  reaction  rate,  but  the  decrease  in 
binding  energy  means  an  increased  desorption  rate,  106. 


percentage  d-band  occupancy 

106 


What  will  be  the  consequence  of  106?  First  of  all  fact  that  the 
groups  are  less  mobile  on  the  surface  of  higher  d  band  filling  may  contribute  to  the 
experimental  fact  that  the  average  hydrocarbon  chain  length  in  the  FT  product 
catalyzed  by  a  metal  at  the  right  side  of  the  transition  series  in  the  Periodic 
Table  is  smaller  than  that  produced  by  a  catalyst  at  the  left  side.5*"*  Secondly  it 
is  experimentally  known  that  the  reactivity  of  an  FT  metal  catalyst  has  a  maximum  as 
one  moves  across  the  transition  series Norskov  attributed  the  trend  to  the 
decrease  in  binding  energy  and  the  adsorption  rate  of  the  adsorbate  along  the 


Table  8.  Binding  Characteristics  of  Ethylene  on  the  Co (0001)  Surface 


Binding 

EteV)5 

Overlap 

Population*3 

1  ' 
Orbital  ‘ 

Occupation  ! 

Geometry 

C-C 

C-M 

a 

TT 

* 

TT 

Total  Charge 
on  Ethylene 

D 

0.83 

0.38 

1.94 

1.75 

1.28 

-1.0 

1 

tt  MimiM  rrttf  **  tt 

1.0 

0.86 

0.23 

1.95 

1.79 

1.08 

1 

-0.4  ' 

- - - 

\  / 

M 

_ 

-0.7 

0.90 

I 

-0.9 

a.  The  binding  energy  is  defined  as  E(slab)  +  E(planar  ethylene)  -  E(system). 
A  positive  sign  means  ethylene  is  bound. 

b.  For  a  comparison  the  C-C  overlap  population  for  the  planar  molecule  is 
1.30  and  for  the  bent  one  1.16. 
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transition  series. ^  The  decrease  in  binding  energy,  he  argued,  is  due  to  an 
increase  in  the  antibonding  between  the  adsorbate  and  the  metal.  This  is  exactly 
what  we  have  in  102.  So  at  the  left  side  of  the  transition  series  the  binding 
is  so  strong  that  it  inhibits  the  coupling  reaction  and  subsequent  desorption.  At 
the  right  side  he  concluded  the  adsorption  rate  is  small  and  the  reactivity  low. 

Our  study  shows  that  the  decrease  in  the  mobility  may  also  contribute  to  the  low 
reactivity  at  the  right  side  of  the  transition  series. 

At  this  point  we  conclude  our  first  look  at  FT  system  in  its  later  stages.  So 
much  more  remains  to  be  done,  and  with  better  calculations  than  ours.  But  we 
believe  that  we  have  gained  some  insight  into  the  essence  of  the  bonding  of  CH^, 

CH2 ,  CH  and  C2H4  to  metal  surfaces,  and  the  migration  propensity  and  coupling 
capabilities  of  these  important  surface  species.  One  useful  conceptual  decom¬ 
position  of  the  barriers,  often  small,  sometimes  large,  that  are  four.c  on  tr.e  *■  iv 
to  products,  is  the  following:  there  are  preferred  sites  of  chemisorption,  dif¬ 
ferential  barriers  to  migration  on  the  surface,  a  proximity  or  crowding  effect 
for  the  nearing  of  fragments  prior  to  reaction,  and,  finally,  an  activation  energy 
for  actual  coupling  and  desorption. 
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APPENDIX 

All  calculations  are  of  the  extended  Huckel  tight  binding  type , ^  with  the 
parameters  given  in  Table  9.  The  M-C  distance  is  chosen  to  be  a  fixed  2 . l8 

throughout  the  calculations  although  experimental  data  indicate  different 

22 

M-C  distances  for  various  adsorption  sites  and  different  C^  fragments.  The  choice 
of  a  constant  M-C  distance  comes  from  our  experience  that  overlap  populations 


for  bonds  of  "unbiased 


or  equal  length  will  be  indicators  of  the  relative  bond 
strengths.  The  Co-Co  distance  in  the  hep  slab  is  2.5l8,^^a  C-H  1.098.  The 
total  energies,  DOS,  COOP  are  calculated  on  a  10-k-point  set  from  ref.  53,  but 
wherever  symmetry  permits  the  10-k-point  set  is  reduced  to  a  7-  or  5-  special  k 


point  set. 


Table  9.  Extended  Huckel  Parameters 


Orbital 

H  .  .  ,  eV 

li 

c2 

a 

C1 

c2 

Co 

3d 

-  9.7 

5.55 

1.9 

0.5448 

0.6556 

4s 

-  7.8 

2.0 

4p 

-  3.8 

2.0 

Cr 

3d 

-  7.9 

4.95 

1.6 

0.4876 

0.7205 

4s 

-  7.3 

1.7 

4p 

-  3.6 

1.7 

Ti 

3d 

-  5.9 

4.55 

1.4 

0.4206 

0.7839 

4s 

-  6.3 

1.5 

4p 

-  3.2 

1.5 

C 

2s 

-21.4 

1.625 

2p 

-11.4 

1.625 

H 

Is 

-13.6 

1.3 

a.  Exponents  and  coefficients  in  a  double  £  expansion  of 
the  3d  orbital. 
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Figure  1:  Total  DOS  (dashed  lines)  of  the  Co(0001) ,  Cr(110)  and  Ti(0001) 

three- layer  slabs.  The  solid  lines  show  the  contribution  from  s  and  p  states  in  a 


and  b  respectively  for  the  Co  case.  indicates  the  Fermi  level.  Cr  and  Ti  total 

state  densities  are  shown  in  c  and  d.  The  d  band  center  of  gravity  and  the  Fermi 
level  shifts  to  higher  energy  on  going  from  Co  to  Ti. 

Figure  2:  A  schematic  picture  showing  the  relative  charge  on  surface  and  bulk  atoms 
(b)  caused  by  different  effective  band  width  for  each  kind  of  atom  (a) . 

Figure  3:  Molecular  orbitals  of  a  CHj  group. 

Figure  4:  Total  DOS  (dashed  line)  and  the  CH^  contribution  (darkened  area)  when  a 
CH^  group  is  chemisorbed  in  an  on-top  geometry  on  Co(0001)  (a)  and  Ti(0001)  (b) 
surfaces.  The  arrows  indicate  the  CH-j  MO  levels  before  the  adsorption  occurs. 

Figure  5:  DOS  of  the  chemisorbed  Co(0001)  (CH-j  on-top),  a.  DOS  before  the 
adsorption,  the  dashed  line  indicates  the  metal  DOS,  the  horizontal  lines  show  the 
free  CH^  MO  levels,  b.  Total  DOS  after  the  chemisorption  occurs.  c.  n-CH^ 
(magnified)  states  in  the  chemisorbed  system.  d  and  e  show  the  dz2  and  s  states 
(magnified)  of  the  metal  atom  below  the  CHj  group. 

Figure  6:  COOP  curve  of  the  M-C  bond  of  a  CH3  group  on  Co(0001)  in  the  on- top 
geometry . 

Figure  7:  A  comparison  of  the  dz2  states,  a  corresponds  to  the  metal  below  the  CH^ 
group  and  b  the  adjacent  metal  atom  (the  one  not  capped  by  a  CH3). 

Figure  8:  COOP  curve  for  the  surface  M-M  bond  of  the  CHj  on  Co(0001)  (on- cop) 
system. 

Figure  9:  Total  DOS  (dashed  line)  and  the  CH-j  contribution  (darkened  area)  of  the 
bridging  CHj  +  Co(0001)  system.  The  arrows  indicate  where  the  free  CH-j  MO  levels 
were  before  the  adsorption. 

Figure  10:  COOP  curves  of  the  two  non-equivalent  C-M  bonds  in  the  bridging  CH-j  + 
Co(0001)  system. 


Figure  11:  DOS  of  the  bridging  CH3  +  Co(0001)  system,  a.  dx2 ;  b.  dz2;  c.  x 
d.  j r  The  dotted  lines  are  the  integrations. 

Figure  12:  DOS  of  the  on-top  CH2  +  Co(0001)  system.  a.  Bare  metal  surface  (dashed 
line)  plus  free  CH2  (MO  levels  indicated  by  horizontal  lines),  b.  Total  DOS  (dashed 
line)  and  the  CH2  contribution  (darkened  area)  for  the  chemisorbed  system.  c.  d,2 
states  of  the  metal  atom  bonded  to  the  CH2  group  (magnified),  d.  s  states 
(magnified),  e.  dyZ  states  (magnified). 

Figure  13:  COOP  curves  of  the  M-M  (solid  line)  and  the  C-M  (dotted  line)  bonds  : 
the  on- top  CH2  +  Co(0001)  system. 

Figure  14:  Total  DOS  (dashed  line)  and  the  CH2  contribution  (darkened  area)  for  the 
bridging  CH2  +  Co(0001)  system  (CH2  perpendicular  to  the  Co-Co  bond). 

Figure  15:  COOP  curves  for  the  M-C  bond.  a.  Perpendicular;  b.  Parallel  CH2  on 
Co(0001) . 

Figure  16:  COOP  curves  for  the  Ml-C  and  M2-C  Bonds.  CH2  group  is  above  the 
triangular  hollow  of  the  Co(0001)  surface. 

Figure  17:  DOS  of  the  CH  +  Co(0001)  slab  system  before  (a)  and  after  (b)  CH 
chemisorbs  on  the  metal  surface  (capping  geometry).  CH  states  are  shown  by 
horizontal  bars  (a)  or  darkened  area  (b) . 

Figure  18:  DOS  of  the  capping  CH  +  Co(0001)  systems,  a.  dz2 ;  b.  dxz  .  Notice 

both  d,2  and  dv_  resonate  with  CH  p  around  -11.2eV. 
c.  az  i  yz 

Figure  19:  n-CH^  DOS  evolution  along  a  coupling  reaction  coordination  9  on 
Co(0001) . 


Figure  20:  C-C  COOP  evolution  along  the  reaction  path  6  on  Co(0001). 

Figure  21:  Relative  energy  of  the  CH3+CH3  system  along  the  coupling  reaction  path 
9.  The  three  curves  correspond  to  different  metal  surfaces  (Co,  Cr  and  Ti) . 

Figure  22:  Methylene  p  state  evolution  along  a  reaction  path  on  Co(0001). 

Figure  23:  Relative  energy  of  the  CH3  +  CH2  system  along  the  reaction  coordinate. 
Three  curves  correspond  to  three  metal  sufaces  (Co,  Cr  and  Ti). 

Figure  24:  Total  DOS  (dashed  line)  and  the  CH2  contribution  (darkened  area)  of  the 
CH2  +  CH2  on  Co(0001)  system. 


IV, 


Figure  25:  A  Walsh  diagram  for  ethylene  from  the  planar  to  the  bent  (30°'  gecm- * r- 
Figure  26:  A  comparison  of  total  DOS  for  ethylene  (bent)  chemisorbed  systems  of 
different  adsorption  geometries  on  Co(0001) .  a)  The  sawhorse  geometry,  b)  The 
butterfly  geometry.  c)  The  capping  geometry.  The  darkened  area  shows  the  ethylene 
contribution. 

Figure  27:  A  comparison  of  the  distribution  of  ethylene  tt*  states  in  three  adsorp¬ 
tion  geometries:  a.  sawhorse;  b.  butterfly;  c.  capping. 


